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R E F A C E 


These  notes  are  based  on  a series  of  lectures  given  in  the 
Lefschetz  Center  for  Dynamical  Systems  in  the  Division  of  Applied 
Mathematics  at  Brown  University  during  the  academic  year  1978-79. 

The  purpose  of  the  lectures  was  to  give  an  introduction  to  the 
applications  of  centre  manifold  theory  to  differential  equations. 

Most  of  the  material  is  presented  in  an  informal  fashion,  by  means 
of  worked  examples  in  the  hope  that  this  clarifies  the  use  of  centre 
manifold  theory. 

The  main  application  of  centre  manifold  theory  given  in  these 
notes  is  to  dynamic  bifurcation  theory.  Dynamic  bifurcation  theory 
is  concerned  with  topological  changes  in  the  nature  of  the  solutions 
of  differential  equations  as  parameters  are  varied.  Such  an  example 
is  the  creation  of  periodic  orbits  from  an  equilibrium  point  as  a 
parameter  crosses  a critical  value.  In  certain  circumstances,  the 
application  of  centre  manifold  theory  reduces  the  dimension  of  the 
system  under  investigation.  In  this  respect  the  centre  manifold 
theory  plays  the  same  role  for  dynamic  problems  as  the  Liapunov- 
Schmitt  procedure  plays  for  the  analysis  of  static  solutions. 

In  order  to  make  these  notes  more  widely  accessible,  we  give 
a full  account  of  centre  manifold  theory  for  finite  dimensional 
systems.  Indeed,  the  first  five  chapters  are  devoted  to  this.  Once 
the  finite  dimensional  case  is  understood,  the  step  up  to  infinite 
dimensional  problems  is  essentially  technical.  Throughout  these  notes 
we  give  the  simplest  such  theory,  for  example  our  equations  are 
autonomous.  Once  the  core  of  an  idea  has  been  understood  in  a simple 


setting,  generalizations  to  more  complicated  situations  are  much  more 
readily  understood. 

In  Chapter  1,  we  state  the  main  results  of  centre  manifold 
theory  for  finite  dimensional  systems  and  we  illustrate  their  use  by 
a few  simple  examples.  In  Chapter  2,  we  prove  the  theorems  which  were 
stated  in  Chapter  1,  and  Chapter  3 contains  further  examples.  In 
Section  2 of  Chapter  3 we  outline  Hopf  bifurcation  theory  for 
2 -dimens ional  systems.  In  Section  3 of  Chapter  3 we  apply  this  theory 
to  a singular  perturbation  problem  which  arises  in  biology.  In 
lixample  3 of  Chapter  6 we  apply  the  same  theory  to  a system  of  partial 
differential  equations.  In  Chapter  4 we  study  a dynamic  bifurcation 
problem  in  the  plane  with  two  parameters.  Some  of  the  results  in  this 
Chapter  are  new  and,  in  particular,  they  confirm  a conjecture  of 
Takens  [49]  . Chapter  4 can  be  read  independently  of  the  rest  of  the 
notes.  In  Chapter  5,  we  apply  the  theory  of  Chapter  4 to  a 4 -dimens ional 
system.  In  Chapter  6,  we  extend  the  centre  manifold  theory  given  in 
Chapter  2 to  a simple  class  of  infinite  dimensional  problems.  Finally, 
we  illustrate  their  use  in  partial  differential  equations  by  means  of 
some  simple  examples. 

I first  became  interested  in  centre  manifold  theory  through 
reading  Dan  Henry's  Lecture  Notes  [30].  My  debt  to  these  notes  is 
enormous.  1 would  like  to  thank  Jack  K.  Hale,  Dan  Henry  and 
John  Mallet-Paret  for  many  valuable  discussions  during  the  gestation 
period  of  these  notes.  In  addition,  my  thanks  go  to  Sandra  Spinacci 
for  converting  my  manuscript  into  this  elegant  typescript. 

This  work  was  done  with  the  financial  support  of  the  United  States 
Army,  Durham,  under  AROD  DAAC  29-76-G0294. 
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Jack  Carr 
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CHAPTER  1 

INTRODUCTION  TO  CENTRE  MANIFOLD  THEORY 


1.1.  Introduct  ion 

In  this  chapter  we  state  the  main  results  of  centre  manifold 
theory  for  finite  dimensional  systems  and  give  some  simple  examples 
to  illustrate  their  application. 


1.2.  Mot i vat  ion 

To  motivate  the  study  of  centre  manifolds  we  first  look  at  a 
simple  example.  Consider  the  system 

k = ax3,  y = -y  + y2,  (1.2.1) 


where  a is  a constant.  Since  the  equations  are  uncoupled  we  can 
easily  show  that  the  zero  solution  of  (1.2.1)  is  asymptotically 
stable  if  and  only  if  a < 0.  Suppose  now  that 


3 2 

x = a x + x y 
y = -y  + y2  xy  - x3. 


(1.2.2) 


Since  the  equations  are  coupled  we  cannot  immediately  decide  if  the 
zero  solution  of  (1.2.2)  is  asymptotically  stable,  but  we  might 
suspect  that  it  is  if  a < 0.  The  key  to  understanding  the  relation 
of  equation  (1.2.2)  to  equation  (1.2.1)  is  an  abstraction  of  the 
idea  of  uncoupled  equations. 
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A curve,  say  y = h(x),  defined  for  |x|  small,  is  said  to 
be  an  invariant  manifold  for  the  system  of  differential  equations 

x * f(x,y) , y = g(x,y)  , (1.2.3) 

if  the  solution  x(t),y(t)  of  (1.2.3)  through  (x0,h(xQ))  lies 
on  the  curve  y = h(x)  for  small  t,  i.e.,  y(t)  = h(x(t)).  Thus, 
for  equation  (1.2.1),  y = 0 is  an  invariant  manifold.  Note  that 
in  deciding  upon  the  stability  of  the  zero  solution  of  (1.2.1),  the 
only  important  equation  is  x = ax  , that  is  we  only  need  study  a 
first  order  equation  on  a particular  invariant  manifold. 

The  theory  that  we  develop  tells  us  that  equation  (1.2.3)  has 

2 

an  invariant  manifold  y = h(x)  , |x|  small,  with  h(x)  = 0(x‘‘)  as 
x -*■  0.  Furthermore,  the  asymptotic  stability  of  the  zero  solution 
of  (1.2.2)  can  be  proved  by  studying  a first  order  equation.  This 
equation  is  given  by 

u = au^  + u^h(u)  = au^  + O(u^),  (1.2.4) 

and  we  see  that  the  zero  solution  of  (1.2.4)  is  asymptotically  stable 
if  a < 0 and  unstable  if  a > 0.  This  tells  us  that  the  zero 
solution  of  (1.2.2)  is  asymptotically  stable  if  a < 0 and  un- 
stable if  a > 0 as  we  expected. 

We  are  also  able  to  use  this  method  to  obtain  estimates  for 
the  rate  of  decay  of  solutions  of  (1.2.2)  in  the  case  a < 0.  For 
example,  if  x(t),y(t)  is  a solution  of  (1.2.2)  with  x(0),y(0) 
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small,  we  prove  that  there  is  a solution  u(t)  of  (1.2.4)  such 
that  x(t)  = u(t) (l  + o(l)) , y (t)  = h(u(t) ) (l+o  (1) ) as  t •+  ®. 

1.3.  Centre  Manifolds 

We  first  recall  the  definition  of  an  invariant  manifold  for 
the  equation 


x ■ N(x)  (1.3.1) 

where  x € lRn.  A set  S c ]Rn  is  said  to  be  a local  invariant 
manifold  for  (1.3.1)  if  for  xQ  € S,  the  solution  x(t)  of  (1.3.1) 
with  x(0)  = xQ  is  in  S for  |t|  <T  where  T > 0.  If  we  can 
always  choose  T = ",  then  we  say  that  S is  an  invariant  manifold. 
Consider  the  system 


x = Ax  + f(x,y) 
9 • By  ♦ g (x , y ) 


(1.3.2) 


where  x € lRn,  y € lRm  and  A and  B are  constant  matrices  such 
that  all  the  eigenvalues  of  A have  zero  real  parts  while  all  the 
eigenvalues  of  B have  negative  real  parts.  The  functions  f and 
g are  C2  with  f(0,0)  - 0,  f * (0,0)  = 0,  g(0,0)  - 0,  g ’ (0,0)  = 0 
( f * is  the  Jacobian  matrix  of  f ) . 

If  f and  g are  identically  zero  then  (1.3.2)  has  two 
obvious  invariant  manifolds,  namely  x = 0 and  y = 0.  The  in- 
variant manifold  x = 0 is  called  the  stable  manifold,  since  if  we 
restrict  initial  data  to  x = 0,  all  solutions  tend  to  zero.  The 
invariant  manifold  y = 0 is  called  the  centre  manifold. 


-4- 


In  general,  if  y = h(x)  is  an  invariant  manifold  for  (1.3.2) 
and  h is  smooth,  then  it  is  called  a centre  manifold  if  h(0)  = 0, 
h'(0)  = 0.  We  use  the  term  centre  manifold  in  place  of  local 
centre  manifold  if  the  meaning  is  clear. 

If  f and  g are  identically  ;ero,  then  all  solutions  of 
(1.3.2)  tend  exponentially  fast,  as  t -*■  ",  to  solutions  of 

x = Ax,  (1.3.3) 

that  is,  the  equation  on  the  centre  manifold  determines  the 
asymptotic  behavior  of  solutions  of  the  full  equation  modulo 
exponentially  decaying  terms.  We  now  give  the  analogue  of  these 
results  when  f and  g are  non -zero.  These  results  are  proved 
in  Chapter  2. 
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' 
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Theorem  1 . There  exists  a centre  manifold  for  (1.3.2),  v 

■> 

| x | < 6 , where  h is  C . 

The  flow  on  the  centre  manifold  is  governed  by  the 
n-d imens ional  system 


h ( x ) , 

m 


u = Au  ♦ f(u,h(u)) 


(1.3.4) 


which  generalizes  the  corresponding  problem  (1.3.3)  for  the  linear 
case.  The  next  theorem  tells  us  that  (1.3.4)  contains  all  the 
necessary  information  needed  to  determine  the  asymptotic  behavior 
of  small  solutions  of  (1.3.2). 


I 


.] 


4 


I 
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Theorem  2 (a)  Suppose  that  the  zero  solution  of  (1.3.4)  is  stable 
(asymptotically  stable)  (unstable).  Then  the  zero  solution  of 
(1.3.2)  is  stable  (asymptotically  stable) (unstable) . 


(h)  Suppose  that  the  zero  solution  of  (1.3.4)  is  stable. 
Let  (x(t),y(t))  be  a solution  of  (1.3.2)  with  (x(0),y(0)) 
sufficiently  small.  Then  there  exists  a solution  u(t)  of  (1.3.4) 
such  that  as  t -*  ®, 


x ( t ) = u(t)  ♦ 0(e  Yt) 
y ( t ) = h(u(t) ) * 0 (e  Yt) 


(1.3.5) 


where  Y > 0 is  a constant. 

If  we  substitute  y(t)  = h(x(t))  into  the  second  equation 
in  (1.3.2)  we  obtain 

h'(x)[Ax  ♦ f ( x , h ( x ) ) ] - Bh(x)  ♦ g(x,h(x)).  (1.3.6) 


liquation  (1.3.6)  together  with  the  conditions  h(0)  = 0,  h'(0)  = 0 
is  the  system  to  be  solved  for  the  centre  manifold.  This  is 
impossible,  in  general,  since  it  is  equivalent  to  solving  (1.3.2). 
The  next  result  however,  shows  that  in  principle,  the  centre  mani- 
fold can  be  approximated  to  any  degree  of  accuracy. 

1 

For  functions  <p : lRn  -*  Ik  which  are  C1  in  a neighborhood 
of  the  origin  define 


(M<P ) (x ) = <p'(x)(Ax  ♦ f(x,4>(x))l  - B4>(x)  - g(x,4>(x)). 
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Note  that  by  (1.3.6),  (Mh)(x)  3 0. 


Theorem  3.  Let  <p  be  a C1  mapping  of  a neighborhood  of  the 
origin  in  lRn  into  lRm  with  $(0)  = 0 and  <f'(0)  = 0.  Suppose 
that  as  x -+  0,  (M4>)(x)  = 0(|x|q)  where  q > 1.  Then  as  x -*■  0 , 
I h C x ) - >J»  (x)  | = 0(  | x|q)  . 

1.4.  Hxarnp les 

We  now  consider  a few  simple  examples  to  illustrate  the  use 
of  the  above  results. 

Hxarnp le  1.  Consider  the  system 

3 , 2 

x = xy  + ax  + by  x 

(1.4.1) 

y = -y  + cx2  ♦ dx“y. 


/ 


* 

a 


i 

Fj 

* -I 


| 


By  Theorem  1,  equation  (1.4.1)  has  a centre  manifold  y « h(x). 
approximate  h we  set 


(M4>)(x)  = ♦,(x)[x*(x)  ♦ ax’*  ♦ bx<P2(x)]  ♦ 4>(x)  " ex1,  - 


dx24>  (x)  . 


- 7- 


If  <P(x)  = 0(x2)  then  (M<f>)(x)  = <P(x)  - cx2  + 0(x4).  Hence,  if 
<P(x)  = cx2,  (M<p)(x)  = 0(x4),  so  by  Theorem  3,  h(x)  = cx2  + 0(x4). 
By  Theorem  2,  the  equation  which  determines  the  stability  of  the 
zero  solution  of  (1.4.1)  is 


I L 


Li 


Q 


u = uh(u)  + au3  + buh2(u)  = (a+c)u3  + 0(u5). 


: 


Thus  the  zero  solution  of  (1.4.1)  is  asymptotically  stable  if 
a + c < 0 and  unstable  if  a + c > 0.  If  a+c=0  then  we 
have  to  obtain  a better  approximation  to  h. 


Suppose  that  a + c = 0.  Let  (x)  = cx2  + v(x)  where 


f(x)  - 0(x4).  Then  (M<P)(x)  = ^(x)  - cdx4  + O(x^).  Thus,  if 


ip(x)  - cx“  + cdx4  then  (Mip)(x)  = O(x^)  so  bv  Theorem  3, 


h(x)  - cx“  + cdx4  + O(x^).  The  equation  that  governs  the  stability 


of  the  zero  solution  of  (1.4.1)  is 


u = uh(u)  + au3  + buh2(u)  = (cd+bc2)u3  + 0(u2) . 


Hence,  if  a + c = 0,  then  the  zero  solution  of  (1.4.1)  is 

2 


asymptotically  stable  if  cd  + be  < 0 and  unstable  if 
2 2 

cd  + be  > 0.  If  cd  + be  = 0 then  we  have  to  obtain  a better 


approximation  to  h (see  Exercise  1). 


L 
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2 

Exercise  1.  Suppose  that  a + c = cd  + be  = 0 in  Example  1. 

Show  that  the  equation  which  governs  the  stability  of  the  zero 
solution  of  (1.4.1)  is  u = •cd‘"u^  + O(u^)  . 


Exercise  2.  Show  that  the  zero  solution  of  (1.2.2)  is  asymptotically 
stable  if  a < 0 and  unstable  if  a > 0. 


Exercise  3.  Suppose  that  in  equation  (1.3.2),  n = 1 so  that  A = 0. 

Suppose  also  that  f(x,y)  = axp  + 0(|x|p+1  ♦ |y|q)  where  2q  > p + 1 
and  a is  non-zero.  Show  that  the  zero  solution  of  (1.3.2)  is 
asymptotically  stable  if  a < 0 and  p is  odd,  and  unstable 
otherwise . 


Example  2.  Consider  the  system 


x = c\  - x*^  ♦ xy 

(1.4.2) 

y = -y  + y"  - X*- 

where  c is  a real  parameter.  The  object  is  to  study  small  solutions 
of  (1.4.2)  for  small  | t | . 
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The  linearized  problem  corresponding  to  (1.4.2)  has  eigenvalues 
-1  and  e.  This  means  that  the  results  given  in  Section  3 do  not 
apply  directly.  However,  we  can  write  (1.4.2)  in  the  equivalent 
form 


• 3 

x = ex  - x + xy 

y = -y  + y2  - x2  (1.4.3) 

e = 0. 

When  considered  as  an  equation  on  ]R3  the  ex  term  in  (1.4.3)  is 
nonlinear.  Thus  the  linearized  problem  corresponding  to  (1.4.3)  has 
eigenvalues  -1,0,0.  The  theory  given  in  Section  3 now  applies  so 
that  by  Theorem  1,  (1.4.3)  has  a two  dimensional  centre  manifold 
y = h(x,e),  |x|  < , | e|  < 6^.  To  find  an  approximation  to  h 

set 


(Mq> ) (x , e)  = <Px(x,  e)  [ex-x3  + x<)>  (x,  e)  ] + <J>(x,e)  + x2  - 4>2(x,e). 

Then,  if  <P  (x)  = -x2,  (M4>)(x,£)  = 0(C(x,£))  where  C is  a homo- 

2 

geneous  cubic  in  x and  e.  By  Theorem  3,  h(x,e)  = -x  + 0(C(x,e)). 
Note  also  that  h(0,e)  = 0 (see  Section  2.6).  By  Theorem  2 the 
equation  which  governs  small  solutions  of  (1.4.3)  is 


u = tu  - 2u3  + 0(|u|C(u,e)) 

e*0. 


(1.4.4, 


The  zero  solution  (u,e)  = (0,0)  of  (1.4.4)  is  stable  for  small  e 


■ 1 ill  winr  - 


so  the  representation  of  solutions  given  by  Theorem  2 applies  here. 

I'or  -6,  < e < 0 the  solution  u = 0 of  the  first  equation  in 
11.4. 4)  is  asymptotically  stable  and  so  by  Theorem  2 the  zero  solution 
of  (1.4.2)  is  asymptotically  stable. 

I'or  0 < e < $2 » solutions  of  the  first  equation  in  (1.4.4) 
consist  of  two  orbits  connecting  the  origin  to  two  small  fixed 
points.  Hence,  for  0 < e < 6?  the  stable  manifold  of  the  origin 
for  (1.4.2)  forms  a separatrix,  the  unstable  manifold  consisting  of 
two  stable  orbits  connecting  the  origin  to  the  fixed  points. 

hxercise  4.  Study  the  behavior  of  all  small  solutions  of 
w + w + ew  + w'^  = 0 for  small  t. 


lixample  3.  Consider  the  equations 

y = -y  + (y  + c)  z 
tz  = y - (y+l)2 


(1.4.5) 


where  e > 0 is  small  and  0 < c < 1.  The  above  equations  arise 
from  a model  of  the  kinetics  of  enzyme  reactions  [29].  If  e = 0, 
then  (1.4.5)  degenerates  into  one  algebraic  equation  and  one 
differential  equation.  Solving  the  algebraic  equation  we  obtain 


z 


(1.4.6) 


and  substituting  this  into  the  first  equation  in  (1.4.5)  leads  to 
the  equation 


where  X * l - c. 

Using  singular  perturbation  techniques,  it  was  shown  in  [29] 
that  for  e sufficiently  small,  under  certain  conditions,  solutions 
of  (1.4.5)  are  close  to  solutions  of  the  degenerate  system  (1.4.6), 

(1.4.7) .  We  shall  show  how  centre  manifolds  can  be  used  to  obtain 
a similar  result. 

Let  t * ei.  We  denote  differentiation  with  respect  to  t by 
• and  differentiation  with  respect  to  i by  '.  Equation  (1.4.5) 
can  be  rewritten  in  the  equivalent  form 

y'  = ef(y,w) 

w'  = - w + y2  - yw  ♦ ef(y,w)  (1.4.8) 
e’  = 0 

where  f(y,w)  » -y  ♦ (y+c)(y-w)  and  w = y - z.  By  Theorem  1, 

(1.4.8)  has  a centre  manifold  w = h(y,c).  To  find  an  approximation 
to  h set 

(M4>)(y,e)  * £4>y  (y , e)  f (y , w)  + h(y,e)  - y2  + yh(y,e)  - ef(y,w). 

If  <Hy,e)  - y2  - Aey  then  (W)(y,c)  * 0(|y|3+|c|3)  so  that  by 
Theorem  3, 


h(y,e)  = y2  - Aey  ♦ 0(|y|3>|e|3). 
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By  Theorem  2,  the  equation  which  determines  the  asymptotic  behavior 
of  small  solutions  of  (1.4.8)  is 

u'  = ef(u,h(u,e)) 

or  in  terms  of  the  original  time  scale 

u * f(u,h(u,e))  * -X(u-u“)  + 0 ( | eu  | ♦ | u | ^ ).  (1.4.9) 

Again,  by  Theorem  2,  if  e is  sufficiently  small  and  y(0),z(0)  are 
sufficiently  small,  then  there  is  a solution  u(t)  of  (1.4.9)  such 
that 


y (t ) = u(t)  ♦ 0(e  t/e) 
z (t)  » y (t ) - h(y(t),e)  ♦ 0(e't/fc). 


(1.4.10) 


Note  that  equation  (1.4.7)  is  an  approximation  to  the  equation  on 

7 

the  centre  manifold.  Also,  from  (1.4.10)  z(t)  - y ( t ) - y (t), 
which  shows  that  (1.4.6)  is  approximately  true. 

The  above  results  are  not  satisfactory  since  we  have  to  assume 
that  the  initial  data  is  small.  In  Chapter  2,  we  show  how  we  can 
deal  with  more  general  initial  data.  Here  we  briefly  indicate  the 
procedure  involved  there.  If  yQ  -1,  then 

(y.w.e)  = (yo.yQd+yo)*1^) 


is  a curve  of  equilibrium  points  for  (1.4.8).  Thus,  we  expect  that 
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there  is  an  invariant  manifold  w *=  h(y,E)  for  (1.4.8)  defined 
for  e small  and  0 < y < m,  (m  = 0(1)),  and  with  h(y,E)  close 


to  the  curve 


2 ,,  . -1 
w * y fl*y) 


(1.4.11) 


For  initial  data  close  to  the  curve  given  by  (1.4.11),  the  stability 
properties  of  (1.4.8)  are  the  same  as  the  stability  properties  of  the 
reduced  equation 

u = f (u,h(u, e) ) . 

1.5.  Bifurcation  Theory 

Consider  the  system  of  ordinary  differential  equations 


w = F (w , e) 
F(0 , e)  = 0 


(1.5.1) 


where  w € 1R  m and  e is  a p-dimensional  parameter.  We  say  that 
e ~ 0 is  a bifurcation  point  for  (1.5.1)  if  the  qualitative  nature 
of  the  flow  changes  at  e = 0,  that  is,  if  in  any  neighborhood  of 
L “ 0 there  exist  points  and  e^  such  that  the  corresponding 

local  phase  portraits  of  (l.S.l)£  and  (l.S.l)e  are  not  topologically 
equivalent . 

Suppose  that  the  linearization  of  (1.5.1)  about  w = 0 is 
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w = C(t)w. 


(1.5.2) 


If  the  eigenvalues  of  C(0)  all  have  non-zero  real  parts  then, 
for  small  |c|,  small  solutions  of  (1.5.1)  behave  like  solutions 
of  (1.5.2)  so  that  e = 0 is  not  a bifurcation  point.  Thus,  from 
the  point  of  view  of  local  bifurcation  theory  the  only  interesting 
situation  is  when  C(0)  has  eigenvalues  with  zero  real  parts. 

Suppose  that  C(0)  has  n eigenvalues  with  zero  real  parts 
and  m eigenvalues  whose  real  parts  are  negative.  We  are  assuming 
that  C(0)  does  not  have  any  positive  eigenvalues  since  we  are 
interested  in  the  bifurcation  of  stable  phenomena. 

Because  of  our  hypothesis  about  the  eigenvalues  of  C(0)  we 
can  rewrite  (1.5.1)  as 


x = Ax  ♦ f (x ,y , e) 

y * By  ♦ g(x,y,c)  (1.5.3) 

t = 0 


where  x 6 lRn , y € lRm,  A is  an  n * n matrix  whose  eigenvalues 
all  have  zero  real  parts,  B is  an  m * m matrix  whose  eigenvalues 
all  have  negative  real  parts,  and  f and  g vanish  together  with 
each  of  their  derivatives  at  (x,y,e)  = (0,0,0). 

By  Theorem  1,  (1.5.3)  has  a centre  manifold  y = h(x,e), 

| x | < 6j,  |e|  < By  Theorem  2 the  behavior  of  small  solutions 

of  (1.5.3)  is  governed  by  the  equation 
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u = Au  ♦ f(u,h(u),e) 
e = 0. 


(1.5.4) 


In  applications  n is  frequently  1 or  2 so  this  is  a very  useful 
reduction.  The  reduction  to  a lower  dimensional  problem  is 
analogous  to  the  use  of  the  Lyapunov- Schmidt  procedure  in  the 
analysis  of  state  problems.  For  the  relationship  between  centre 
manifold  theory  and  other  perturbation  techniques  such  as  amplitude 
expansions,  see  [12]. 

1.6.  Comments  on  the  Literature 

Theorems  1-3  are  the  simplest  such  results  in  centre  manifold 
theory  and  we  briefly  mention  some  of  the  possible  generalizations. 

(1)  The  assumption  that  the  eigenvalues  of  the  linearized 
problem  all  have  non-positive  real  parts  is  not  necessary. 

(2)  The  equations  need  not  be  autonomous. 

(3)  In  certain  circumstances  we  can  replace  'equilibrium 
point'  by  'invariant  set'. 

(4)  Similar  results  can  be  obtained  for  certain  classes  of 
infinite-dimensional  evolution  equations,  such  as  partial  differ- 
ential equations. 

There  is  a vast  literature  on  invariant  manifold  theory 
[1,6,18,19,23,24,26,28,30,31,36,37,38,40,42].  For  applications  of 
invariant  manifold  theory  to  bifurcation  theory  see  [1,12,13,14,15, 
20,27,30,32,33,34,39,40,41,42,45,50].  For  a simple  discussion  of 
stable  and  unstable  manifolds  see  [18,  Chapter  13]  or  [23,  Chapter  3]. 


In  Chapter  2 we  prove  Theorems  1-3.  Our  proofs  of  Theorems  1 
and  2 are  modelled  on  Kelly  [37,38].  Theorem  3 is  a special  case  of 
a result  of  Henry  [30]  and  our  proof  follows  his.  Throughout 
Chapter  2 we  use  methods  that  generalize  to  infinite  dimensional 
problems  in  an  obvious  way. 


CHAPTER  2 

PROOFS  OF  THEOREMS 


2.1.  Introduction 

In  this  chapter  we  give  proofs  of  the  three  main  theorems 
stated  in  Chapter  1.  The  proofs  are  essentially  applications  of 
the  contraction  mapping  principle.  The  procedure  used  for  defining 
the  mappings  is  rather  involved,  so  we  first  give  a simple  example 
to  help  clarify  the  technique.  The  proofs  that  we  give  can  easily 
be  extended  to  the  corresponding  infinite  dimensional  case;  indeed 
essentially  all  we  have  to  do  is  to  replace  the  norm  |*  | in  finite 
dimensional  space  by  the  norm  ||‘  ||  in  a Banach  space. 

2.2.  A Simple  Example 

We  consider  a simple  example  to  illustrate  the  method  that  we 
use  to  prove  the  existence  of  centre  manifolds. 

Consider  the  system 

*1  = x2,  x,  = 0,  >v  « -y  ♦ g(x1,x2),  (2.2. U 

where  g is  smooth  and  g(x  ,x,)  - O(x^x^)  as  (x  .xj  - (0,0). 

We  prove  that  (2.2.1)  has  a local  centre  manifold. 

2 

Let  4/:  1R~  -►  1R  be  a C00  function  with  compact  support  such 
that  ^(Xj ,x,)  a 1 for  (x  ,x^)  in  a neighborhood  of  the  origin. 
Define  G by  Glx^.x,)  * f (x ^ , x 1) g (x^ , x ) . We  prove  that  the 
system  of  equations 
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x^  = Xj,  x2  ■ 0,  y = -y  + G(xx>x2),  (2.2.2) 

2 

has  a centre  manifold  y = h(x^,x2),  (x^.x^)  £ IR  . Since 

G(xx,x7)  = g(Xj,x2)  in  a neighborhood  of  the  origin,  this  proves 

2 2 

that  y = h(x^,x2),  x^  + x2  < 6 for  some  6,  is  a local  centre 
manifold  for  (2.2.1). 

The  solution  of  the  first  two  equations  in  (2.2.2)  is 
xL  (t)  = Zj  + z7t,  x2(t)  = z2,  where  Xj,(0)  = z^.  If 

y(t)  = h(Xj (t)  ,x7(t) ) is  a solution  of  the  third  equation  in  (2.2.2) 
then 


jy  h(zx+z2t,z2)  = -h(z1+z2t,z2)  + G(z1+z2t,z2).  (2.2.3) 

To  determine  a centre  manifold  for  (2;2.2)  we  must  single  out  a 
special  solution  of  (2.2.3).  Since  G(xx,x2)  is  small  for  all  x^ 
and  x7,  solutions  of  the  third  equation  in  (2.2.2)  behave  like 
solutions  of  the  linearized  equation  y = -y.  The  general  solution 
of  (2.2.2)  therefore  contains  a term  like  e t.  As  t -*■  °°,  this 
component  approaches  the  origin  perpendicular  to  the  z\>z2  Plane- 
Since  the  centre  manifold  is  tangent  to  the  z\*zi  plane  at  the 
origin  we  must  eliminate  the  e t component,  that  is  we  must 
eliminate  the  component  that  approaches  the  origin  along  the  stable 
manifold  as  t -*  00 . To  do  this  we  solve  (2.2.3)  together  with  the 
condit  ion 


lim  h(zx+z2t , z2)et 
£ ->  - 00 


= 0. 


(2.2.4) 


Integrating  (2.2.3)  between 


- oo 


and  0 and  using  (2.2.4) 


we  obtain 


r° 

h(zltz2)  = j esC.(z1  + z7s,z2)ds. 

By  construction,  y = h(Zj,z7)  is  an  invariant  manifold  for  (2.2.2). 
Using  the  fact  that  G has  compact  support  and  that  G(x^,x,)  has 
a second  order  zero  at  the  origin  it  follows  that  h(Zj,z?)  has  a 
second  order  zero  at  the  origin, that  is  h is  a centre  manifold. 


2.3.  Ex i stence  of  Centre  Manifolds 

In  this  section  we  prove  that  the  system 


* = Ax  ♦ f ( x , y ) 
X * By  ♦ g (x , y ) 


(2.3.1) 


has  a centre  manifold.  As  before  x £lRn,  y € lRm,  the  eigenvalues 
of  A have  zero  real  parts,  the  eigenvalues  of  B have  negative 
real  parts  and  f and  g are  C*"  functions  which  vanish  together 
with  their  derivatives  at  the  origin. 

Theorem  1 . liquation  (2.3.1)  has  a local  centre  manifold 
y = h(x),  | x | < 6,  where  h is  C". 

Lroo f . As  in  the  example  given  in  the  previous  section,  we  prove 
the  existence  of  a centre  manifold  for  a modified  equation.  Let 
y:  lRn  •+  (0,1)  be  a C°°  function  with  Y(x)  = 1 when  |x|  < 1 
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and  V(x)  = 0 when  |x|  > 2.  For  e > 0 define  F and  G by 
F(x,y)  = f (x<K|)  ,y) , G(x,y)  = g(xV(|),y). 


The  reason  that  the  cut-off  function  y is  only  a function  of  x 
is  that  the  proof  of  the  existence  of  a centre  manifold  generalizes 
in  an  obvious  way  to  infinite  dimensional  problems. 

We  prove  that  the  system 


x = Ax  ♦ F(x,y) 
y = By  ♦ G(x,y) 


(2.3.2) 


has  a centre  manifold  y = h(x),  x € lRn,  for  small  enough  e.  Since 
F and  G agree  with  f and  g in  a neighborhood  of  the  origin, 
this  proves  the  existence  of  a local  centre  manifold  for  (2.3.1). 

For  p > 0 let  X be  the  set  of  Lipschitz  functions 
h : lRn  -»  lRm  with  Lipschitz  constant  p,  |h(x)|  < p for  x £ 
and  h(0)  = 0.  With  the  supremum  norm  ||'||,  X is  a complete 
space . 

For  hex  and  x(1  € lRn , let  x(t,Xg,h)  be  the  solution  of 


x = Ax  + F (x , h ( x) ) , x (0 , xQ , h)  = xQ.  (2.3.3) 


The  bounds  on  F and  h ensure  that  the  solution  of  (2.3.3)  exists 
for  all  t.  We  now  define  a new  function  Th  by 


(Th)(x0)  = 


-Bs 


G(x(s,x0,h),h(x(s,x0,h)))ds.  (2.5.4) 
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If  h is  a fixed  point  of  (2.3.4)  then  h is  a centre  manifold 
for  (2.3.2).  We  prove  that  for  p and  e small  enough,  T is  a 
contraction  on  X. 

Using  the  definitions  of  F and  G,  there  is  a continuous 
function  k(e)  with  k ( 0 ) = 0 such  that 

I F (x  ,y) | + | G(x,y) | < k(c), 

I F (x , y)  - F (x ' ,y ' ) | < k(e)[|x-x'|  + |y-y'|],  (2.3.5) 

I G (x , y ) - G (x ' , y ' ) | < k( e) [ | x-x ' | + | y-y  * | ] , 

for  all  x,  x’  € lRn  and  all  y,  y’  € lRm  with  |y|,  |y'|  < e 
Since  the  eigenvalues  of  B all  have  negative  real  parts, 
there  exist  positive  constants  3,C  such  that  for  s < 0 and 
y € lRm , 


|e'Bsy|  < Ce 6s | y | . (2.3.6) 

Since  the  eigenvalues  of  A all  have  zero  real  parts,  for  each 
r > 0 there  is  a constant  M(r)  such  that  for  x € lRn  and 
s € 1R, 

| eAsx | < M(r)er'sl |x| . (2.3.7) 

Note  that  in  general,  M(r)  -*•  00  as  r -*•  0 . 

If  p < £,  then  we  can  use  (2.3.5)  to  estimate  terms  involving 
G(x(s,x0,h) ,h(s,xQ,h) ) and  similar  terms.  We  shall  suppose  that 
p < £ from  now  on. 
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If  xn  € 1R  , then  using  (2.3.6)  and  the  estimates  on  G and 
h,  we  have  from  (2.3.4)  that 


| Th(xn) | < C6'1k(e). 


(2.3.8) 


Now  let  Xq,  x^  £ 1R Using  (2.3.7)  and  the  estimates  on  F 
and  h,  we  have  from  (2.3.3)  that  for  r > 0 and  t < 0, 


x ( t , x Q , h ) - x(t ,Xj ,h) | < M(r)e'rt|xQ-x1 

r(s-t) 


+ 2M(r)k(e) 


e 1 J|x(s,xQ,h)  - x (s  f x j , h) |ds 


By  Gronwall's  inequality,  for  t < 0, 


(t,xQ,h)  - x (t , Xj  , h)  | < M(r)  Ixj-Xgle'^, 


(2.3.9) 


where  Y - r + 2M(r)k(t).  Using  (2.3.9)  and  the  bounds  on  G and 
h,  we  obtain  from  (2.3.4) 


|Th(xQ)  - Th ( x x ) | < C ( M ( r ) + p) k ( e) ( 0- y) ' 1 | xQ -x} | (2.3.10) 


if  e.  and  r are  small  enough  so  that  0 > y. 

Similarly,  if  hj,  h2  € X and  xQ  €lRn,  we  obtain 


l lhl (xQ)  • Th2(x0)|  < Ck(e)r'1(0-y)'1[4M(r)k2(e)  + 11  l|h1-h2|| . (2.3.11) 


. mmiM  Mia 


.... - 
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By  a suitable  choice  of  p,e  and  r,  we  see  from  (2.3.8), 

(2.3.10)  and  (2.3.11)  that  T is  a contraction  on  X.  This  proves 

the  existence  of  a Lipschitz  centre  manifold  for  (2.3.2).  To  prove 

that  h is  C*  we  show  that  T is  a contraction  on  a subset  of  X 

consisting  of  Lipschitz  differentiable  functions.  The  details  are 

similar  to  the  proof  given  above  so  we  omit  the  details.  To  prove 
2 

that  h is  C we  imitate  the  proof  of  Theorem  4.2  on  page  333 

of  [18] . 

2.4.  Reduction  Principle 

The  flow  on  the  centre  manifold  is  governed  by  the  n-dimensional 
system 


u - Au  ♦ f(u,h(u)).  (2.4.1) 

In  this  section  we  prove  a theorem  which  enables  us  to  relate  the 
asymptotic  behavior  of  small  solutions  of  (2.3.1)  to  solutions  of 
(2.4.1). 

We  first  prove  a lemma  which  describes  the  stability  properties 
of  the  centre  manifold. 

Lemma  1 . Let  (x(t),y(t))  be  a solution  of  (2.3.2)  with  | (x(0) ,y (0) ) | 
sufficiently  small.  Then  there  exist  positive  and  p such  that 

I y ( t ) - h(x(t) ) | < Cle'vt  | y (0)  - h(x(0))| 


for  all  t > 0. 
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Proof  . Let  (x(t),y(t))  be  a solution  of  (2.3.2)  with  (x(0),y(0)) 
sufficiently  small.  Let  z(t)  = y(t)  - h(x(t)),  then  by  an  easy 
computation 


z = Bz  + N(x,z) 


(2.4.2) 


where 


N(x,z)  = h' (x) [ F (x , h (x) ) - F (x , z + h (x) ) ] ♦ G (x , z + h (x) ) - G(x,h(x)). 


Using  the  definitions 

of  F 

and 

G and  the 

bounds  on  h, 

there  is 

a continuous  function 

6(e) 

with 

6(0)  = 0 

such  that 

1 N (x , z) | < 6(c) | z | if 

1 z 1 

< e. 

Using  (2.3. 

6)  we  obtain, 

from 

(2.4.2), 

I z ( t ) | < Ce" et | z (0) | + C6(e)  [ V 3 (t ‘ s) | z (s ) | ds 

1 0 

and  the  result  follows  from  Gronwall's  inequality. 

Before  giving  the  main  result  in  this  section  we  make  some 
remarks  about  the  matrix  A.  Since  the  eigenvalues  of  A all  have 
zero  real  parts,  by  a change  of  basis  we  can  put  A is  the  form 
A = Aj  + A7  where  A2  is  nilpotent  and 

At 

I e x|  = | x | . (2.4.3) 

Since  A^  is  nilpotent,  we  can  choose  the  basis  such  that 


! 


I 

i 


1 

1 

T 

i 


1. 


p 


1 
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I A2x|  < (6/4) | x | , (2.4.4) 

where  0 is  defined  by  (2.3.6). 

We  assume  for  the  rest  of  this  section  that  a basis  has  been 
chosen  so  that  (2.4.3)  and  (2.4.4)  hold. 

Theorem  2.  (a)  Suppose  that  the  zero  solution  of  (2.4.1)  is  stable 

(asymptotically  stable) (unstable) . Then  the  zero  solution  of  (2.3.1) 
is  stable  (asymptotically  stable) (unstable)  . 

(b)  Suppose  that  the  zero  solution  of  (2.3.1)  is  stable. 
Let  (x(t),y(t))  be  a solution  of  (2.3.1)  with  (x(0),y(0)) 
sufficiently  small.  Then  there  exists  a solution  u(t)  of  (2.4.1) 
such  that  as  t -*• 


1. 


r 

L 


x ( t ) = u ( t ) ♦ 0(e*Yt) 
y(t)  = h(u(t) ) ♦ 0 (e  Yt ) 


(2.4.5) 


[ 
• « 

m 

L 

l 

I 

I 

I 


where  Y > 0 is  a constant  depending  only  on  B. 

Proof.  If  the  zero  solution  of  (2.4.1)  is  unstable  then  by 
invariance,  the  zero  solution  of  (2.3.1)  is  unstable.  From  now  on 
we  assume  that  the  zero  solution  of  (2.3.1)  is  stable.  We  prove 
that  (2.4.5)  holds  where  (x(t),y(t))  is  a solution  of  (2.3.2)  with 
| (x  (0)  ,y(0))  | sufficiently  small.  Since  F and  G are  equal  to  f 
and  g in  a neighborhood  of  the  origin  this  proves  Theorem  2.  We 
divide  the  proof  into  two  steps. 


i 

I 


I 


I 

I- 


‘ 


' 


j 

' S 


! 
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I.  Let  Uq  £ ]Rn  and  Zg  € lRm  with  |(Ug,Zg)|  sufficiently 
small.  Let  u(t)  be  the  solution  of  (2.4.1)  with  u(0)  = Ug.  We 
prove  that  there  exists  a solution  (x(t),y(t))  of  (2.3.2)  with 
y(0)  - h(x(0))  = Zg  and  x(t)  - u(t),  y(t)  - h(u(t))  exponentially 
small  as  t -*•  °°. 

II.  By  Step  I we  can  define  a mapping  S from  a neighborhood 
of  the  origin  in  lRn  + m into  ]Rn*m  by  S(Ug,Zg)  = (Xg,Zg)  where 
Xg  = x(0).  For  | (Xg,Zg)|  sufficiently  small,  we  prove  that 
(Xg,Zg)  is  in  the  range  of  S. 

I.  Let  (x(t),y(t))  be  a solution  of  (2.3.2)  and  u(t)  a 
solution  of  (2.4.1).  Note  that  if  u(0)  is  sufficiently  small, 

u = Au  ♦ F (u , h (u) ) (2.4.6) 

since  solutions  of  (2.4.1)  are  stable.  Let  z(t)  = y(t)  - h(x(t)), 
'P(t)  = x(t)  - u(t),  then  by  an  easy  computation 

z * Bz  + N (4>  +u , z)  (2.4.7) 

= A<1>  + R(4>,z)  (2.4.8) 

where  N is  defined  in  the  proof  of  Lemma  1 and 

R (4>  , z ) = F(u*4>  , z + h(u+4> ) ) - F(u,h(u)). 


We  now  formulate  (2.4.7),  (2.4.8)  as  a fixed  point  problem. 


For  a > 0,  K > 0,  let  X be  the  set  of  continuous  functions 


V : 10,-)  -*■  lRn  with  |4>(t)eat|  < K for  all  t > 0.  If  we  define 
| | 4>  | | = sup{  | (t)  eat  | : t > 0},  then  X is  a complete  space.  Let 
u^.Zq  be  sufficiently  small  and  let  u(t)  be  the  solution  of 
(2.4.0)  with  u(0)  = uQ.  Given  £ X let  z(t)  be  the  solution 
of  (2.4.7)  with  z(0)  = z^.  Define  T<}>  by 


r°°  A.  (t-s) 

(T4> ) ( t ) = - j e 1 (A2<J>(s)  ♦ R(<P(s)  ,z(s))]ds.  (2.4.9) 


l. 


1. 


L 


i. 


r* 


*• 


We  solve  (2.4.9)  by  means  of  the  contraction  mapping  principle. 

If  ♦ is  a fixed  point  of  T,  then  retracing  our  steps  we  find  that 
x ( t ) = u ( t ) + <P(t),  y(t)  = z ( t ) + h(x(t) ) is  a solution  of  (2.3.2). 
We  can  take  a to  be  as  close  to  3 as  we  please  at  the  cost  of 
increasing  K and  shrinking  the  neighborhood  on  which  the  result 
is  valid.  For  simplicity  however,  we  take  K = 1 and  2a  = 3 
where  3 is  defined  by  (2.3.6). 

Using  the  bounds  on  F,G,h  and  the  fact  that  N(4>,0)  = 0,  there 
is  a continuous  function  k(c)  with  k(0)  = 0 such  that  if 
G and  zi>z7  € ]Rm  with  | z ^ | < e,  then 


IN^j.Zj)  - N(*2,z2)|  < k(e)  ( | z2  | -♦2  | ♦ | Zj  - z2  | J 

l. 

^ - R(^>2,22)|  5 k(c)[|zrz2|  ♦ | 4>  x 2 | 1 . 

I 

I 
I 


» 


Pf 


' I 


From  (2.4.7) , 
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2(t)|  5 C I z0  I e'  Bt  * Ck(e)  |V  0(t ' s)  | z (s)  | ds 


where  we  have  used  (2.3.6)  and  (2.4.10).  By  Gronwall's  inequality 


-B.t 

z(t)|  < C|z|e  1 


(2.4.11) 


where  = 3 - Ck(£).  From  (2.4.9),  if  c is  sufficiently 


sma 1 1 , 


-B,s 


|'R(t)  | < | ♦ k(e)  j (e  as  ♦ C|z0|e  1 


) ds  < e 


where  we  have  used  (2.4.3),  (2.4.4),  (2.4.10)  and  (2.4.11).  Hence 
T maps  X into  X. 

Now  let  e * and  let  z\,z2  0e  ^'e  coresP°nding 

solutions  of  (2.4.4)  with  z^O)  = :0.  We  first  estimate 
w(t)  = z . ( t ) - z,(t).  From  (2.4.7)  , 


'(t)|  < Ck  ( t)  [te'0(t's)[|z,(s)||4.1(s)  - q>  ,(s 
n 11 


s) I ♦ |w(s)  Ids. 


Using  (2.4.11), 


w(t)|  < Cjk(t)  | |e‘0t  ♦ Ck(t)  j e'  0(t's1  |w(s)  |ds 


where  is  a constant,  so  that  by  Gronwall's  inequality 


lwlt)|  = |Zj(t)  - 22  ( t ) I < Cjk(e)  I |$1*4>:|  |e  1 . (2.4.12) 


-29- 


Using  (2.4.4)  and  (2.4.12),  for  £ sufficiently  small, 

|Wj(t)  - T*2(t)|  < | ll^-^U  + k(e)  ^(^(s)  - 02(s)  | + | Zj  (s)  - z2(s)|)ds 

< o|K1-^2|| 


where  u < l . 

The  above  analysis  proves  that  for  each  (uQ,z0)  sufficiently 
small  , T has  a unique  fixed  point.  If  U is  a neighborhood  of 
the  origin  in  lRn+m  then  it  is  easy  to  repeat  the  above  analysis 
to  show  that  T : X * U -*•  X is  a continuous  uniform  contraction. 

This  proves  that  the  fixed  point  depends  continuously  on  u^ 
and  Zq. 

II.  Define  S by  S(u0,zQ)  = (x0,zQ)  where  xQ  = uQ  + <*> ( 0 ) . 
Since  depends  continuously  on  uQ  and  zQ,  S is  continuous. 

We  prove  that  S is  one-to-one,  so  that  by  the  Invariance  of 
Domain  Theorem  (see  [9]  or  (47])  S is  an  open  mapping.  Since 
S(0 , 0)  = 0,  this  proves  that  the  range  of  S is  a full  neighbor- 
hood of  the  origin  in  lRn  + m. 

Proving  that  S is  one-to-one  is  clearly  equivalent  to  proving 
that  if  uQ  + <PQ  C 0 ) = Uj  + 4>  L ( 0 ) then  uQ  = Uj  and  <PQ(0)  = >^(0). 
If  uQ  + 4>q(0)  = u^  + 4^(0)  then  the  initial  values  for  x and 
y are  the  same,  so  that  by  uniqueness  of  solution  of  (2.3.2), 
uQ(t)  ♦ <pQ(t)  = u^(t)  + <j>  j ( t ) for  all  t > 0,  where  u^(t)  is  the 
solution  of  (2.4.6)  with  u^(0)  = u^.  Hence,  for  t > 0, 

u0(t)  - Uj(t)  = $x(t)  - 4>0(t). 


(2.4.13) 


Since  the  real  parts  of  the  eigenvalues  of  A are  all  zero, 
lim  |u.(t)  - un(t)|eLt  - ® for  any  c > 0 unless  u. (0)  = un(0). 

t-*oo  1 U 1 II 

Also,  |vp  j(t)  | < e at  for  all  t > 0.  It  now  follows  from  (2.4.13) 
that  S is  one-to-one  and  this  completes  the  proof  of  the  theorem. 

2 . 5 . Approximation  of  the  Centre  Manifold 

for  functions  : lRn  ■*  lRm  which  are  in  a neighborhood  of 

the  origin  define 

(M'P)U)  = 4>'(x)[Ax  ♦ f(x,<J>(x)))  - B4>(x)  - g(x,<J>(x)). 

Theorem  3.  Suppose  that  4>(0)  = 0,  4> ' (0)  = 0 and  that  (M4>)(x)  = 
0(|x|q)  as  x - 0 where  q > 1.  Then  as  x + 0, 

I h (x)  - Mx)|  = 0(|x|q). 

Proo_f.  Let  d;  lRn  -*  lRm  be  a continuously  differentiable  function 
with  compact  support  such  that  0 (x)  = ^(x)  for  |x|  small.  Set 

N(x)  = 0 ' (x) [Ax  ♦ F(x ,0 (x) ) ] - B0(x)  - G(x,0(x)),  (2.S.1) 

where  1:  and  G are  defined  in  Theorem  1.  Note  that  N(x)  = 

0(|x|q)  as  x -*■  0 . 

In  Theorem  1,  we  proved  that  h was  the  fixed  point  of  a 
contraction  mapping  T:  X -*  X.  Define  a mapping  S by  Sz  = 

ng  a closed  subset  Y c X,  Since 
, S is  a contraction  mapping  on  Y. 


' 
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For  K > 0 let 

Y = {z  € X:  j z (x)  | < K|x|q  for  all  x £]Rn}. 

If  we  can  find  a K such  that  S maps  Y into  Y then  we  will 
have  proved  the  theorem. 

We  first  find  an  alternative  formulation  of  the  map  S.  For 
z £ Y let  x(t,Xg)  be  the  solution  of 


x = Ax  + F(x,z(x)  + 0(x)),  x(0,xn)  = x (0 ) . (2.5.2) 


From  (2.3.4) 


(T  ( z+ti  ) ) (xn)  = 


e'BsG(x(s,x0) , z(x(s,xQ))  + 0 (x(s ,xQ) ) )ds 


Now 


-e(x0)  - - 


[e'Bs0(x(s,xo))]ds 


e'Bs[B0(x(s,xo))  - ^ 0(x(s,xo))Jds. 


Writing  x for  x(s,xQ)  etc.,  from  (2.5.1)  and  (2.5.2) 


B°(x)  - 3j0(x)  = Bti(x)  - 0 ' (x)  [Ax  + F (x , z (x)  + 0(x))] 


-N(x)  - G(x,0(x))  + 0 ' (x) (F(x,0)  - F(x,z(x)  + d(x))] 


. 
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Using  Sz  = T(z+9)  - 9 and  the  above  calculations 

r°  -bs 

(Sz)(xQ)  * j Q(x(s,x0),  z(x(s,x0)))ds  (2.5.3) 

where  x(s,x0)  is  the  solution  of  (2.5.2)  and 

Q(x,z)  * G(x,9*z)  - G(x,9)  - N(x)  ♦ 9'(x)(F(x,9)  - F(x,9+Z))  . (2.5.4) 

We  now  show  that  S maps  Y into  Y for  some  K > 0.  By 
choosing  o suitably,  we  may  assume  that  |Q(x)|  < e for  all 
x € JRn.  Since  N(x)  = 0(|x|q)  as  x - 0, 

I N(x)  | < Cj  | x | q , x € lRn  (2.5.5) 

where  is  a constant.  Now 

| Q (x , z ) | < | Q(x , 0) | + |Q(x,z)  - Q (x , 0) | 

(2.5.6) 

= I N (x) | ♦ | Q(x , z)  - Q(x , 0) | . 


We  can  estimate  |Q(x,z)  - Q(x,0)|  in  terms  of  the  Lipschitz 
constants  of  F and  G.  Using  (2.3.5),  there  is  a continuous 
function  k(£)  with  k(0)  * 0,  such  that 

lQ(x,z)  - Q (x , 0) | < k(e) | z | (2.5.7) 

for  | z|  < e.  using  (2.5.S),  (2.5.6),  (2.5.7),  for  z € Y and 
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and  x 6 lRn,  we  have  that 


I Q (x , z) | < Cx | x | q 4 k(c)|z(x)| 

< (Cj+Kk (fc) ) ) x | q . 


(2.5.8) 


Using  the  same  calculations  as  in  the  proof  of  Theorem  1,  if 
x(t,xQ)  is  the  solution  of  (2.5.2),  then  for  each  r > 0,  there  is 
a constant  M(r)  such  that 


|x(t,x0)|  < M(r) |x0|e'Yt,  t < 0 (2.5.9) 

where  Y = r ♦ 2M(r)k(t). 

Using  (2.3.6),  (2.5.8)  and  (2.5.9),  if  z £ Y, 

I (Sz)  (xQ)  ( < C(C1>kk(c))(M(r))t<(B-qY)'1|x0|q  = C2  | xQ  | q 

provided  e and  r are  small  enough  so  that  3 - qY  > 0.  By  choosing 

K large  enough  and  e small  enough,  we  have  that  C-,  < K and  this 

completes  the  proof  of  the  theorem. 

2.6.  Properties  of  Centre  Manifolds 

(1)  In  general  (2.3.1)  does  not  have  a unique  centre  manifold, 

lor  example,  the  system  x = -x\  y = -y,  ha*  the  two  parameter 

family  of  centre  manifolds  y = hU.Cj.c.,)  where 
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h (x , c ,c,)  = 


1 - 2 

c.  exp(-  T x ) , 


l*c2)  = | 0 


, 1 -2, 

2 exP('  J X 5 


x > 0 

x = 0 
x < 0 


However,  if  h and  are  two  centre  manifolds  for  (2.3.1), 

then  by  Theorem  3,  h(x)  - h^(x)  = 0(|x|^)  as  x -►  0 for  all 

q > 1. 

(2)  If  f and  g are  C^,  (k  > 2),  then  h is  [37]. 

If  f and  g are  analytic,  then  in  general  (2.3.1)  does  not  have 
an  analytic  centre  manifold,  for  example  it  is  easy  to  show  that 
the  system 

‘ 3 • 2 

x = -x  , y * -y  ♦ x (2.6.1) 

does  not  have  an  analytic  centre  manifold  (see  exercise  (1)). 

(3)  Centre  manifolds  need  not  be  unique  but  there  are  some 

points  which  must  always  be  on  any  centre  manifold.  For  example, 
suppose  that  is  a smaH  equilibrium  point  of  (2.3.1)  and 

let  y = h(x)  be  any  centre  manifold  for  (2.3.1).  Then  by  Lemma  1 
we  must  have  y^  = h(x^).  Similarly,  if  T is  a small  periodic 
orbit  of  (2.3.1),  then  T must  lie  on  all  centre  manifolds. 

(4)  Suppose  that  x(t),y(t)  is  a solution  of  (2.3.1)  which 
remains  in  a neighborhood  of  the  origin  for  all  t > 0.  An 
examination  of  the  proof  of  Theorem  2,  shows  that  there  is  a solution 


■I 


r 
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u(t)  of  (2.4.1)  such  that  the  representation  (2.4.5)  holds. 

(5)  In  many  problems  the  initial  data  is  not  arbitrary,  for 
example,  some  of  the  components  might  always  be  nonnegative. 

Suppose  S c IRn+m  with  0 € S and  that  (2.3.1)  defines  a local 
dynamical  system  on  S.  It  is  easy  to  check,  that  with  the  obvious 
modifications,  Theorem  2 is  valid  when  (2.3.1)  is  studied  on  S. 

Exercise  1.  Consider 

• 3 • 2 

x=-x,  y = -y  + x . (2.6.1) 

Suppose  that  (2.6.1)  has  a centre  manifold  y = h(x) , where  h is 
analytic  at  x = 0.  Then 

h(x)  . lax" 
n=  2 

for  small  x.  Show  that  a,,  = 0 for  all  n and  that  a = na 

2n+l  n+2  n 

for  n = 2,4,...,  with  a,  = 1.  Deduce  that  (2.6.1)  does  not  have 
an  analytic  centre  manifold. 

Exercise  2 (Modification  of  an  example  due  to  S.J.  van  Strien  [48]). 

If  f and  g are  C functions,  then  for  each  r,  (2.3.1) 
has  a C centre  manifold.  However,  the  size  of  the  neighborhood 
on  which  the  centre  manifold  is  defined  depends  on  r.  The  following 
example  shows  that  in  general  (2.3.1)  does  not  have  a C°°  centre 
manifold,  even  if  f and  g arc  analytic. 
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Cons  ider 


x 


- tx  - 


-y 


-*• 


e 


0. 


(2.6.2) 


oo 

Suppose  that  (2.6.2)  has  a C centre  manifold 
| x | < <5,  1 1 1 < 6.  Choose  n > 6 *.  Then  since 

oo 

C in  x,  there  exist  constants  a^,a, a1n 


y = h(x,c)  for 
h(x,(2n)  * ) is 
such  that 


h ( x , ( 2 n ) '*) 


rn  i 2n+lx 

l a.x  + 0(x  ) 

i=  1 1 


for  | x | small  enough.  Show  that  a^  = 0 for  odd  i and  that  if 
n > 1 , 


(1  - (2i)(2n)'1)a2.  = (2i-2)a,._2>  i = 2,...,n 


(2.6.3) 


a,  ^ 0. 


Obtain  a contradiction  from  (2.6.3)  and  deduce  that  (2.6.2)  does 

oo 

not  have  a C centre  manifold. 


Exercise  3.  Suppose  that  the  nonlinearities  in  (2.3.1)  are  odd, 
that  is  f (x , y ) = -f(-x,-y),  g(x,y)  = -g(-x,-y).  Prove  that  (2.3.1) 
has  a centre  manifold  y = h(x)  with  h(x)  = -h(-x).  (The  example 
x = -x\  y = -y,  shows  that  if  h is  any  centre  manifold  for  (2.3.1) 
then  h(x)  / -h(-x)  in  general.] 


2.7.  Global  Invariant  Manifolds  for  Singular  Perturbation  Problems 


To  motivate  the  results  in  this  section  we  reconsider 


Example  3 in  Chapter  1.  In  that  example  we  applied  centre  manifold 


theory  to  a system  of  the  form 


y'  ■ ef(y*w) 


w'  = - w + y“  - yw  + tf(y,w) 


(2.7.1) 


£ ' =0 


where  f(0,0)  = 0.  Because  of  the  local  nature  of  our  results  on 


centre  manifolds,  we  only  obtained  a result  concerning  small  initial 


data.  Let  v = -w(l+y)  + y , then  we  obtain  a system  of  the  form 


y'  = ^jCy.v) 


V = -V  ♦ eg2(y,v) 


(2.7.2) 


£ ' = 0 


where  g^(0,0)  = 0,  i = 1,2.  Note  that  if  y / -1,  then  (y,0,0) 
is  always  an  equilibrium  point  for  (2.7.2)  so  we  expect  that  (2.7.2) 
has  an  invariant  manifold  v = h(y,e)  defined  for  -1  < y < m say 


and  t.  sufficiently  small. 


Theorem  4 . Consider  the  system 


x ' = Ax  ♦ ef (x ,y , t) 


y'  = By  ♦ tg(x,y,£) 


(2.7.3) 


£ ' « 0 
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where  x £ IRn , y € IRm  and  A,B  are  as  in  Theorem  1.  Suppose  also 
that  f,g  are  C"  with  f(0,0,0)  = 0,  g (0,0,0)  = 0.  Let  m > 0. 
Then  there  is  a 6 > 0 such  that  (2.7.3)  has  an  invariant  manifold 
y = h(x,0,  | x | < m,  |t|  < 6,  with  |h(x,e)|  < C|c|,  where  C is  a 
constant  which  depends  on  m,A,B,f  and  g. 

Proof . Let  y : lRn  -*  {0,11  be  a C function  with  y(x)  = 1 if 
| x | < m and  y(x)  = 0 if  |x|  > m ♦ 1.  Define  F and  G hy 

I'lx.y.c)  = ef(xy(x) ,y,t) , G(x,y,t)  = »-g  (xY  (x)  ,y  , c)  . 

We  can  then  prove  that  the  system 


x'  = Ax  ♦ F(x,y,t) 
y'  = By  ♦ G(x,y,c) 


.4) 


has  an  invariant  manifold  y = h(x,t),  x € lRn,  for  |t|  sufficiently 
small.  The  proof  is  essentially  the  same  as  that  given  in  the  proof 
of  Theorem  1 so  we  omit  the  details. 


Kema  r k . If  x = (Xj.x-, x(i)  then  we  can  similarly  prove  the 

existence  of  h(x,t)  for  m < x.  < nt.. 

1 ’ J -l  i i 

The  flow  on  the  invariant  manifold  is  given  hy  t he  equation 


u ' = Au  ♦ tf  t.u,h(u,c))  . 


(2.7.3) 
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With  the  obvious  modifications  it  is  easy  to  show  that  the 
stability  of  solutions  of  (2.7.3)  is  determined  by  equation  (2.7.5) 
and  that  the  representation  of  solutions  given  in  (2.4.5)  holds. 
Finally,  we  state  an  approximation  result. 

Theorem  5.  Let  <p : ]Rn+1  ]Rm  satisfy  $(0,0)  = 0 and 
|(M$)(x,e)|  < CeP  for  |x[  < m where  p is  a positive  integer, 

C is  a constant  and 

(M$)(x,£)  « Dx$(x,£)[Ax  + e f(x,4>(x,e))J  - B$(x,c)  - e g(x  ,$  (x,  e) ) . 
Then,  for  |x|  < m. 


I h(x , e)  - 4»(x,e)|  < C 


for  some  constant  C^. 

Theorem  5 is  proved  in  exactly  the  same  way  as  Theorem  3 so  we 
omit  the  proof. 

For  further  information  on  the  application  of  centre  manifold 
theory  to  singular  perturbation  problems  see  Fenichel  120]  and 
Henry  [ 30] . 


CHAPTER  3 
EXAMPLES 


3.1.  Rate  of  Decay  Estimates  in  Critical  Cases 

In  this  section  we  study  the  decay  to  zero  of  solutions  of 
the  equation 


r ♦ r ♦ f (r)  = 0 (3.1.1) 

where  f is  a smooth  function  with 

f(r)  = r3  + ar’  ♦ 0(r7)  as  r -*■  0 , (3.1.2) 

where  a is  a constant.  By  using  a suitable  Liapunov  function  it 
is  easy  to  show  that  the  zero  solution  of  (3.1.1)  is  asymptotically 
stable.  However,  because  f'(0)  = 0,  the  rate  of  decay  cannot  be 
determined  by  linearization. 

In  (8)  the  rate  of  decay  of  solutions  was  given  using 
techniques  which  were  special  to  second  order  equations.  We  show 
how  centre  manifolds  can  be  used  to  obtain  similar  results. 

We  first  put  (3.1.1)  into  canonical  form.  Let  x = r + r, 
v = r,  then 


x = -f(x-y) 
y = -y  - f (x-y) . 


(3.1.3) 


By  Theorem  1 of  Chapter 


1 

* 9 


(3.1.3)  has  a centre  manifold  v = h(x) . 
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By  Theorem  2 of  Chapter  2,  the  equation  which  determines  the 
asymptotic  behavior  of  small  solutions  of  (3.1.3)  is 


u = -f(u-h(u)) 


(3.1.4) 


Using  (3.1.2)  and  h(u)  = 0(u4)I 


u = -u3  + 0(u4) 


(3.1.5) 


Without  loss  of  generality  we  can  suppose  that  the  solution  u(t) 
of  (3.1.5)  is  positive  for  all  t > 0.  Using  [.'Hospital's  rule, 


-1  = Lim^y=  Lim  t’ 1 
t-wo  u'’  t-*-°° 


ru(t) 


s 3ds . 


Hence,  if  w(t)  is  the  solution  of 


w = -w  , w(0)  = 1, 


(3.1.6) 


then  u ( t ) = w(t+o(t)).  Since 


w(t)  = — t"1^"  + Ct-3/2  + nr*-'5/2' 


/I 


0 ( t *"fc) 


(3.1.7) 


where  C is  a constant,  we  have  that 


u(t)  = — t'1/2  ♦ o(t‘1/2). 

/I 


(3.1.8) 


To  obtain  further  terms  in  the  asymptotic  expansion  of  u(t), 


we 
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need  an  approximation  to  h(u) . To  do  this,  set 

(M<j> ) (x)  = -<J>  ’ (x)  f (x-4>  (x) ) + )>(x)  * f(x-4>(x)). 

If  4>(x)  = -x3  then  (J#)(x)  = 0(x3)  so  that  by  Theorem  3 
of  Chapter  2,  h(x)  = -x3  + 0(x3).  Subsituting  this  into  (3.1.4) 
we  obtain 


u = -u3  - (a+3)u3  + 0(u2)  . 


(3.1.9) 


Choose  T so  that  u(T)  = 1.  Dividing  (3.1.9)  by  u3,  integrating 
over  [T,t]  and  using  (3.1.8),  we  obtain 


w 1(u(t))  = t + constant  + (3  + a) 


f\.2 


u (s)ds  (3.1.10) 


where  w is  the  solution  of  (3.1.6).  Using  (3.1.8)  and  (3.1.9), 


u (s)ds  = 
J T 


^ sT  + 

\j  1 

-1/2 

= -In  t + constant  + o(l). 


(3.1.11) 


Substituting  (3.1.11)  into  (3.1.10)  and  using  (3.1.7), 


u ( t ) = — t'1/2  - 1 


-3/2 


n 


4/2 


( (a+  3)  In  t + C]  + o(t'3/2)  (5.1.12) 


where  C is  a constant. 

If  x(t),y(t)  is  a solution  of  (3.1.3),  it  follows  from 
Theorem  2 of  Chapter  2 that  either  x(t),y(t)  tend  to  zero  ex- 
ponentially fast  or 
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x ( t ) = ±u(t) , y(t)  = +u3(t) 
where  u(t)  is  given  by  (3.1.12). 

3.2.  Hopf  Bifurcation 

There  is  an  extensive  literature  on  Hopf  Bifurcation  [1,13,15, 

16,27,30,35,39,40,42,46],  so  we  only  given  an  outline  of  the  theory. 

Our  treatment  is  based  on  [15] . 

Consider  the  one-parameter  family  of  ordinary  differential 
_ 2 

equations  on  1R  , 


x = f (x,a)  , x € 1R2, 

such  that  f(0,a)  = 0 for  all  sufficiently  small  u.  Assume  that 
the  linearized  equation  about  z = 0 has  eigenvalues  y(u)  ± iio(u) 
where  Y(0)  = 0,  w(0)  = w * o.  We  also  assume  that  the  eigenvalues 
cross  the  imaginary  axis  with  nonzero  speed  so  that  y'(0)  * 0. 

Since  Y'(0)  * 0,  by  the  implicit  function  theorem  we  can  assume 
without  loss  of  generality  that  Y C01 ) = u*  By  means  of  a change  of 
basis  the  differential  equation  takes  the  form 

* = A(ot)x  + F(x,^)  , (3.2.1) 


where 


A(a) 


a -w(u) 
ui(a)  a 


F(x,u)  = 0( | x | 2)  . 
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I 
I 

Under  the  above  conditions,  there  are  periodic  solutions  of  j| 

(3.2.1)  bifurcating  from  the  zero  solution.  More  precisely,  for 
small  there  exists  a unique  one  parameter  family  of  small  amplitude  | 

periodic  solutions  of  (3.2.1)  in  exactly  one  of  the  cases 

(i)  u < o,  (ii)  u = o,  (iii)  u > 0.  However,  further  conditions  * 

on  the  nonlinear  terms  are  required  to  determine  the  specific 

i 

type  of  bifurcation. 

Hxercise  1.  Use  polar  co-ordinates  for 

x ^ - wx,  + KXj (Xj+x“) 

2 2 

X,  = wXj  + + Kx^(xJ+xp 

to  show  that  case  (i)  applies  if  K > 0 and  case  (iii)  applies  if 
k < 0. 

To  find  periodic  solutions  of  (3.2.1)  we  make  the  substitution 
Xj  = er  cos  x = tr  sin  f u -♦  ta)  (3.2.2) 


where  is  a function  of  u.  After  substituting  (3.2.2)  into  (3.2.1) 

we  obtain  a system  of  the  form 


r = t(ar  ♦ r”CT  (d  ,at)  ] + *-“r^C.(0  ^*t)  O(t^) 

4 

J + 0(0. 


(3.2.3) 


I 


J 


We  now  look  for  periodic  solutions  of  (3.2.3)  with  t.  -*  0 and  r 
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near  a constant  r^.  If  Cj  and  are  independent  of  0 and 
the  higher  order  terms  are  zero  then  the  first  equation  in  (5.2.5) 
takes  the  form 


i • £[«r  ♦ gr2]  + e2Kr3  . (5.2.4) 

Periodic  solutions  are  then  the  circles  r = r^,  where  r^  is  a 
zero  of  the  right  hand  side  of  (5.2.4).  We  reduce  the  first 
equation  in  (5.2.5)  to  the  form  (5.2.4)  modulo  higher  order  terms 
by  means  of  a certain  transformation.  It  turns  out  that  the 
constant  6 is  zero.  Under  the  hypothesis  K is  non-zero,  it  is 
straightforward  to  prove  the  existence  of  periodic  solutions  by 
means  of  the  implicit  function  theorem.  The  specific  type  of 
bifurcation  depends  on  the  sign  of  K so  it  is  necessary  to  obtain  a 
formula  for  K. 

Let  F(x,u)  = lFj(x1,x2,a)  , F, (Xj ,x  , , a)J T and  let 

F‘j(xl,x2,a)  = BJ(x1,x2,u)  + fii  (Xj  ,x2 , '•*)  + 0(xJ  + x1)  (5.2.5) 

where  B-j  is  a homogeneous  polynomial  of  degree  i in  (x^.x,). 
Substituting  (5.2.2)  into  (5.2.1)  and  using  (3.2.5)  we  obtain 
(3.2.3)  where  for  i = 3,4, 

C^0.01)  = (cos  b ) B j _ j (cos  0,  sin  u,a) 

(sin  b)B2_j(cos  b,  sin  <3,u). 


♦ 


(3.2.0) 


-46- 


Lemma  1.  There  exists  a coordinate  change 


r = r + tUj  (r  ,6  ,u  , e)  + e u,  (r  ,0  ,u , e) 


which  transforms  (3.2.3)  into  the  system 


r = tar  + t“r^K  + 0(e3) 


0 = + 0(t) 


where  the  constant  K is  given  by 


K = (l/2i>) 


2" 


lc4(w,o)  - woic,(y ,o)D3(6,o))de 


where  Cj  and  are  given  by  (3.2.6) 


and 


(3.2.7) 


(3.2.8) 


D,(o ,0)  = cos  0B,(cos  0,  sin  0,0)  - (sin  0)B*(cos  0,  sin  0,0) 


The  coordinate  change  is  constructed  via  averaging.  We  refer 
to  I 15]  for  a proof  of  the  Lemma. 

If  K = 0 then  we  must  make  further  coordinate  changes.  We 
assume  that  K * 0 from  now  on. 

Recall  that  we  are  looking  for  periodic  solutions  of  (3.2.7) 
with  t -*  0 and  r near  a constant  r^.  This  suggests  that  we  set 
a = -sgn(K)t  and  r(1  = |K|  The  next  result  gives  the 

existence  of  periodic  solutions  of 


r = t“[-sgn(K)r  ♦ r3k)  ♦ 0(t3) 

o = u>0  + 0(t) 


(3.2.9) 


i 
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with  r - Tq  small. 

Lemma  2.  liquation  (3.2.9)  has  a unique  periodic  solution  for  £ 
small  and  r in  a compact  region  either  for  £ > 0 (when  K < 0) 
or  t.  < 0 (when  K > 0)  . Also 

r = rQ  + 0(£),  <3(t,£)  = u»Qt  + 0 ( £ ) 

and  the  period  of  the  solution  t(e)  is  given  by 

T(£)  = (2ir/u>0)  ♦ 0 ( £)  . 

The  periodic  solution  is  stable  if  K < 0 and  unstable  if  K > 0. 

Lemma  2 is  proved  by  a simple  application  of  the  implicit  func- 
tion theorem.  We  again  refer  to  (15]  for  a proof. 

Lemma  2 also  proves  the  existence  of  a one  parameter  family  of 
periodic  solutions  of  (3.2.1).  We  cannot  immediately  assert  that 
this  family  is  unique  however,  since  we  may  have  lost  some  periodic 
solutions  by  the  choice  of  scaling,  i.e.  by  scaling  m + at  and 
by  choosing  £ = -sgn(K)ot. 

In  order  to  justify  the  scaling  suppose  that  = R cos  0, 

x,  = R sin  0 is  a periodic  solution  of  (3.2.1)  bifurcating  from 
Xj  = x,  = 0.  Then  R satisfies 

R = aR  + 0(R2) . 

When  R attains  its  maximum,  R = 0 so  that  R = 0(a).  This 


. * 


-48- 


justifies  the  scaling  u -*•  ue.  A similar  argument  applied  to 
periodic  solutions  of  (3.2.7)  justifies  the  choice  of  e. 

Theorem . Suppose  that  the  constant  K defined  by  (3.2.8)  is  non- 
zero. Then  (3.2.1)  has  a unique  periodic  solution  bifurcating  from 
the  origin,  either  for  u > 0 (when  K < 0)  or  u < 0 (when 
K > 0) . If  x = R cos  0,  y = R sin  then  the  periodic  solution 
has  the  form 


R(t,<*)  = l-K-1!172  + 0 ( | ot  | ) 

«(t,u)  = wQt  ♦ 0(|a|1/2) 

with  period  t(u)  = (21'/Ujq)  + 0(|u|^^2).  The  periodic  solution  is 
stable  if  K < 0 and  unstable  if  K > 0. 

Finally,  we  note  that  since  the  value  of  K depends  only  on 
the  nonlinear  terms  evaluated  at  u = 0 , when  applying  the  above 
Theorem  to  (3.2.1)  we  only  need  assume  that  the  eigenvalues  of  Av.a) 
cross  the  imaginary  axis  with  non-zero  speed  and  that 


with  Wq  non-zero. 

3.3.  Hopf  Bifurcation  in  a Singular  Perturbation  Problem 

In  this  section  we  study  a singular  perturbation  problem  which 
arises  from  a mathematical  model  of  the  immune  response  to 


i i 


> 


<: 


ij 


11 

f' 
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ant igen  143).  The  equations  arc 

lx  * -lx3  ♦ (a  - ^)x  ♦ b - \ ] 
a * 6(1 -x)  - a - Vjab  (3.3.1) 

b = - Y . a b ♦ > b 

I *, 

whore  1-,6,y.,Y,  are  positive  parameters.  In  the  above  model  a 
and  b represent  certain  concentrations  so  they  must  be  non 
negative.  Also,  x measures  the  stimulation  of  the  system  and  it 
is  scaled  so  that  |x|  s 1.  The  stimulation  is  assumed  to  take 
place  on  a much  faster  time  scale  than  the  response  so  that  < is 
very  sma 1 1 . 

The  above  problem  was  studied  in  143)  and  we  briefly  outline  t he 
method  used  by  Merrill  to  prove  the  existence  of  periodic  solutions 
of  (3.3.1).  Putting  l = o in  the  first  equation  in  (3.3.1)  we 
oht a i n 

x3  ♦ (a  - y)x  ♦ b - ^ * 0 • (3.3.2) 


Solving  (3.3.2)  for  x as  a function  of  a and  h we  obtain  x 1 (a.h) 
and  substituting  this  into  the  second  equation  in  (3.3.1)  we  obtain 


a ■ v 6 ( 1 I (a , b ) ) a > j ah 
b - -Yjftb  ♦ Y ,h. 


(3.3.3) 


■^samsmmtaa 


I 
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Using  <5  as  the  parameter,  it  was  shown  that  relative  to  a certain 
equilibrium  point  ( a ^ , h ^ ) a Hop f bifurcation  takes  place  in 
(3.3.3).  Bv  appealing  to  a result  in  singular  perturbation  theory, 
it  was  concluded  that  for  e sufficiently  small,  (3.3.1)  also  has 
a periodic  solution. 

We  use  the  theory  given  in  Chapter  2 to  obtain  a similar  result 
Let  (Xy.apjb^)  be  a fixed  point  of  (3.3.1).  If  b()  * 0 then 

ao  = and  va0  satisfy 


X 


3 

0 


1 

T 


V?b 


2 0 


0. 


(3.3.4) 


Recall  also  that  for  the  biological  problem,  we  must  have  b(1  > 0 

and  |x(J  <1.  We  assume  for  the  moment  that  x(1  and  b()  satisfy 

(3.3.4)  and  these  restrictions.  The  reality  of  these  solutions  are 

considered  later.  We  let  a ^ for  the  rest  of  this  section 

Let  y = a - a(),  z ■ b - b(),  w = -y(x-x(1)  - x()v  - 2 where 
,2  1 

v = ^xo  + »0  ' T • Then  assuming  y is  non- zero, 


c\V  = g(w,y,z,c) 

y = f,(w,y,z,e)  (3.3.5) 

i = f 3 ( w , y , z , t ) 


where 
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g(w,y,z,e)  = fj(wfy,z,e)  - £x0f2(w,y,z,£)  - £f3(w,y,z,£) 
fjCw.y.z.e)  = -vw  ♦ N (w+xQy+z ,y) 

f 2 (w»  y » 2 » e)  = (|  ^'lxo  - 1 - Yibo)y  + ^ 1 ' Y2)z  + § ^"lw  " Y!yz 

fj(w,y,z,e)  = * Yjb0y  - Yjyz 

N(t),y)  = -4>’2b3  + 3^’lx002  - y6. 

In  order  to  apply  centre  manifold  theory  we  change  the  time 
scale  by  setting  t = £s.  We  denote  differentiation  with  respect  to 
s by  ’ and  differentiation  with  respect  to  t by  • . Equation 

(3.3.5)  can  now  be  written  in  the  form 

W = g(wfy,z,e) 
y’  = tf2(w, y,z,e) 

(3.3.6) 

z*  = ef3(w,y,z,e) 

£'  = 0. 


Suppose  that  v > 0.  Then  the  linearized  system  corresponding  to 

(3.3.6)  has  one  negative  eigenvalue  and  three  zero  eigenvalues.  By 
Theorem  1 of  Chapter  2,  (3.3.6)  has  a centre  manifold  w = h(y,zft-). 
By  Theorem  2 of  Chapter  2,  the  local  behavior  of  solutions  of  (3.3.6) 
is  determined  by  the  equation 


y'  * £f2(h(y,z,£) ,y,z,t) 
Z'  = ef3(H(y,z,£) ,y,z,t) 


(3.3.7) 


or  in  terms  of  the  original  time  scale 
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y = f2(h(y,z,e) ,y,z,c) 
z = f3(h(y,z,t) , y , z , t) . 


(3.3.8) 


We  now  apply  the  theory  given  in  the  previous  section  to  show  that 
(3.3.8)  has  a periodic  solution  bifurcating  from  the  origin  for 
certain  values  of  the  parameters. 

The  linearization  of  the  vector  field  in  (3.3.8)  about 
y = z = 0 is  given  by 


J(e) 


Ylb0 


*Ylb0 


+ 0(C). 


If  (3.3.8)  is  to  have  a llopf  bifurcation  then  we  must  have 
trace(J(c))  = o and  y v - y~  > 0.  from  the  previous  analysis, 
we  must  also  have  that  x^ , b()  are  solutions  of  (3.3.4)  with 
|x()|  < 1,  b„  > 0 and  > 0.  We  do  not  attempt  to  obtain  the 
general  conditions  under  which  the  above  conditions  are  satisfied, 
we  only  work  out  a special  case. 

Lemma . Let  Yj  < 2 y , . Then  for  each  e > 0,  there  exists  6(c), 

x()(<0,  b0(t)  such  that  0 < 2x0U)  < 1,  bQU)  > 0,  * > 0, 
cS(c)v  * - 2y7  > 0,  trace  J(t)  = 0 and  (3.3.4) 


is  satisfied. 
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Proof.  Fix  and  Y,  with  Y^  < 2Y2*  If  xO’*50,<5  satisfy  the 

second  equation  in  (3.3.4)  then 

y 

trace ( J (0) ) = j (xQv'1  - ^ (l-x0)l 

It  is  easy  to  show  that  there  is  a unique  xn(0)  € (0,-^-)  that 
satisfies  trace(J(0))  = 0.  Clearly  ^ > 0 for  this  choice  of 
x(1(0).  We  now  obtain  b{1(0)  and  6(0)  as  the  unique  solution  of 
(3.3.4)  and  an  easy  computation  shows  that  bg(0)  > 0,  6(0)  > 0 
and  6(0)Y  * - 2y2  > 0.  By  the  implicit  function  theorem,  for 
c,  Xg  - x^(O),  - bg(0)  sufficiently  small,  there  exists 

6(c,xt1,b0)  = 6(0)  + 0(c)  such  that 

trace  (J (c) ) = | ' 1 xQ  - 1 - Yjbg  + 0(c)  = 0. 


After  substituting  6 = 6(e,x^,b^)  into  (3.3.4),  another  applica- 
tion of  the  implicit  function  theorem  gives  the  result.  This  completes 
the  proof  of  the  Lemma. 

From  now  on  we  fix  Yj  and  y2  with  Y j < 2y~,.  Using  the 
same  calculations  as  in  the  Lemma,  for  each  t and  6 with  t 
and  6 - 6(c)  sufficiently  small  there  is  a solution  x(1(c,<5), 
b0(t,6)  of  (3.3.4).  Writing  x^  = xQ(e,6(c))  and  trace(J)  as  a 
function  of  6,  we  have  that 


|y  (trace (J  16))) 


6=6(l) 


irry-  (<,YlX°  * 


(2Y2-Y1)x() 


l 2 Y j x 0 ) ♦ 0(c)  < 0, 
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if  £ is  sufficiently  small.  Hence,  the  eigenvalues  of  J(6) 
cross  the  imaginary  axis  with  non-zero  speed  at  6 = 6(e). 

Let  y j = fc(e)z,  z = m(e)y  where 

H£)  = (Y1bQ(e))'1/2  ♦ 0(e),  m(e)  = [ ( 6/2)>P ' 1 - Y2)1/2  + 0(e). 

For  6 = 6(e),  (3.3.8)  in  these  new  coordinates  becomes 

h ’ “Vi  - Y1m"1C£)y1zi 

S1  = “oyl  * m(E)(5/2)*,'1h(m'1(E)z1>  «'1(E)y1,e)  - Yj11' 1 (EJyjZj 
where 


«J(0  = Ytb0(e)  [(6/2)V_1  - Y2)  ♦ 0(e)  . 

To  apply  the  results  of  the  last  section,  we  need  to  calculate 
the  K(e)  associated  with  (3.3.9).  We  shall  show  how  K(0)  can  be 
calculated;  if  K(0)  is  non-zero  then  K(t)  will  be  non-zero  also. 
To  calculate  K(0)  we  need  to  know  the  quadratic  and  cubic  terms  in 
(3.3.9)  when  e = o.  Thus,  we  have  to  find  h(y,z,0)  modulo  fourth 
order  terms. 

Let 


(M9)(y,z)  = -g(<P  (y,  z)  ,y,  z,0)  (3.3.10) 


Then  by  Theorem  2 of  Chapter  2,  if  we  can  find  <P  such  that 
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4.  4 


(M<P)(y,z)  = 0(y*,  + zH)  then  h(y,z,0)  = 4>(y,z)  ♦ 0(y4+z4).  Suppose 


that 


(p  +4 

*2  ^ 3 


(3.3.11) 


where  is  a homogeneous  polynomial  of  degree  j.  Substituting 

(3.3.11)  into  (3.3.10)  we  obtain 


(M*)(y.z)  = V4>2  (y , z)  - ^'^(Xgy+z)2  + y(xQy+z)  + 0(|y|3  ♦ |z|3) 


Hence,  if 


(P2(y.z)  3 3^  " ^Xq (xQy+  z)  2 - V ' Xy (xQy+z) 


(3.3.12) 


then  (M4»)(y,z)  = 0(|y|3  ♦ |z|3).  Substituting  (3.3.11)  into 
(3.3.10)  with  2 (y , z)  given  by  (3.3.12),  we  obtain 


(M$> ) (y , z)  = V<J>3(y,z)  ♦ v'2(X()y+z)3-  6V_1x  Q(x  y+z)*  2(y,z) 


4.  4, 


♦ y4,2(y,z)  ♦ 0 (y  *z  ) . 


Hence , 


h(y,z)  = *2(y,z)  - v ' J(x0y+z)3+  6Y'2x0(x0yz)*2(y,  z) 


* V 1 y^  2 ( y , Z ) ♦ 0(y4  + z4) 


1 1 u 

s 

1 


where  <f’2(y,z)  is  defined  by  (3.3.12).  K(0)  can  now  be  calculated 
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as  in  the  previous  section.  The  sign  of  K(0)  will  depend  on  y} 
and  Y , • If  K(0)  is  non-zero  then  we  can  apply  Theorem  1 of 
Section  2 to  prove  the  existence  of  periodic  solutions  of  (3.3.1). 
The  stability  of  the  periodic  solutions  is  determined  by  the  sign  of 
K'(0).  If  K(0)  is  zero  then  we  have  to  calculate  K(e). 


CHAPTER  4 

BIFURCATIONS  WITH  TWO  PARAMETERS  IN  TWO  DIMENSIONS 


4.1.  Introduction 

In  this  chapter  we  consider  an  autonomous  ordinary  differ- 
ential equation  in  the  plane  depending  on  a two-dimensional 
parameter  c.  We  suppose  that  the  origin  x = 0 is  a fixed  point 
for  all  c.  More  precisely,  we  consider 


x = f(x,e),  x € IR2,  e = (e1,e2)  eiR2, 
f(0,e)  = 0. 


(4.1.1) 


The  linearized  equation  about  x = 0 is 

x = A(e)x, 


and  we  suppose  that  A(0)  has  two  zero  eigenvalues.  The  object  is 

to  study  small  solutions  of  (4.1.1)  for  (£.,£-)  in  a full 

neighborhood  of  the  origin.  More  specifically,  we  wish  to  divide 

a neighborhood  of  t-  = 0 into  distinct  components,  such  that  if 

are  in  the  same  component,  then  the  phase  portraits  of  (4.1.1) 

and  (4.1.1) are  topologically  equivalent.  We  also  want  to  describe 

c 

the  behavior  of  solutions  for  each  component.  The  boundaries  of 
the  components  correspond  to  bifurcation  points. 

Since  the  eigenvalues  of  A(0)  are  both  zero  we  have  that 
either  (i)  A(0)  is  the  zero  matrix,  or  (ii)  A(0)  has  a Jordan  block, 


A ( 0) 


x 


- 


There  is  a distinction  between  (i);nd  (ii)  even  for  the  study  of 
fixed  points.  Under  generic  assumptions,  in  Case  (ii),  equation 
(4.1.1)  has  exactly  2 fixed  points  in  a neighborhood  of  the  origin. 
Tor  Case  (i)  the  situation  is  much  more  complicated  (17,  25]. 

Another  distinction  arises  when  we  consider  the  eigenvalues  of 
A(e).  We  would  expect  the  nature  of  the  eigenvalues  of  A(c")  to 
determine  (in  part)  the  possible  types  of  bifurcation.  If 


A(t)  = 


£1  0 


then  the  eigenvalues  of  A(t-)  are  always  real  so  we  do  not  expect 
to  obtain  periodic  orbits  surrounding  the  origin.  On  the  other 
hand  i f 


(4.1.2) 

then  the  range  of  the  eigenvalues  of  A(t)  is  a neighborhood  of 
the  origin  in  I,  that  is,  if  ; is  a small  complex  number  then 
A(t)  has  an  eigenvalue  : for  some  t. 

We  shall  assume  from  now  on  that  A(e)  is  given  hv  (4. 1.2.). 


A(t-)  = 


0 1 


£1  L2 


■ 


* 
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Takens  (49,  5 0 J and  Bodganov  [see  2]  have  studied  normal 
forms  for  local  singularities.  Takens  shows,  for  example,  that  any 
perturbation  of  the  equation 

*1  = x2  + xi * x2  = X1 
is  topologically  equivalent  to 

X1  = x2  + xl’  x2  = el  + e2xl  ' X1 

for  some  There  arc  certain  difficulties  in  applying  these 

results  since  we  must  transform  our  equation  into  normal  form,  modulo 
higher  order  terms. 

In  [39,  p.  333-348],  Kopell  and  Howard  study  (4.1.1)  under  the 
assumptions  that  A(e)  is  given  by  (4.1.2)  and  that 


° f •> 


where  f 1 is  the  second  component  of  f.  Their  approach  consists 
of  a systematic  use  of  scaling  and  applications  of  the  implicit 
function  theorem. 

In  this  chapter,  we  use  the  same  techniques  as  Kopell  and 
Howard  to  study  (4.1.1)  when  the  nonl inearitics  arc  cubic.  Our 
results  confirm  the  conjecture  made  by  Takens  [49]  on  the  bifurcation 
set  of  (4.1.1). 

The  results  on  quadratic  nonlinearities  are  given  in  Section  9 
in  the  form  of  exercises.  Most  of  these  results  can  be  found  in 
Kopell  and  Howard  [39]. 

4.2.  Prel iminar ies 

Consider  equation  (4.1.1)  where  A(t)  is  given  by  (4.1.2).  We 
also  suppose  that  the  linearization  of  f(x,t)  is  A(e)x,  and  that 
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f(0,t)  = 0,  f(x,c)  = -£(*x,e). 


(4.2.1) 


The  object  is  to  study  the  behavior  of  all  small  solutions  of  (4.1.1) 

for  t in  a full  neighborhood  of  the  origin.  Equation  (4.1.1)  is 

still  too  general,  however,  so  we  shall  make  some  additional  hv- 

T 

potheses  on  the  nonlinear  terms.  Set  f = (f^.f^)  , 


3x; 


f,(o,o),  e 


r 


T f,(o,m 


axjx,' 


(HI) 

u * 

0 

(H2) 

6 * 

0 

(H3) 

35 

fl 

(0,0)  = 0. 

3x  j 

(HI) 

imp  lies 

that  for  small 

t ixed  points.  Under  (HI),  it  is  easy  to  show  that  by  a change  of 
co-ordinates  in  (4.1.1),  we  can  assume  (H5)  (see  Remark  1). 

We  assume  (H3)  in  order  to  simplify  the  computations. 

Under  (HI)  we  can  prove  the  existence  of  families  of  periodic 
orbits  and  homoclinic  orbits.  Under  (H1)-(H3)  we  can  say  how  many 
periodic  orbits  of  (4.1.1)  exist  for  fixed  t.  The  sign  of  8 will 
determine  the  direction  of  bifurcation  and  the  stability  of  the 
periodic  orbits  among  other  things. 

from  now  on  we  assume  (H1)-(H3). 

The  main  results  arc  given  in  Figures  2-5.  Sections  3-8  of 
this  chapter  show  how  we  obtain  these  pictures.  The  pictures  for 

8 > 0 are  obtained  by  using  the  change  of  variables  x,  * -x,, 

, * -t*,  and  t -*  -t. 


REGION  I 


REGION  2 


REGION  3 


REGION  4 


-ot>- 


The  cases  u > 0 and  ot  < 0 are  geometrically  different. 

The  techniques  involved  in  each  case  are  the  same  and  we  only  do  the 
more  difficult  case  u < o . The  case  u > o is  left  for  the 
reader  as  an  exercise. 

From  now  on  we  assume  u < 0 in  addition  to  (HI) -(M3).  Note 
that  this  implies  that  locally,  (4.1.1)  has  1 fixed  point  for 
< 0 and  3 fixed  points  for  > 0. 

4.3.  Scaling 

We  scale  the  variables  in  equation  (4.1.1)  so  that  the  first 
components  of  the  non-zero  fixed  points  are  given  by  ±1  + 0(t). 

To  do  this  we  introduce  parameters  u,6,  scaled  variables  y^.y., 
and  a new  time  t by  the  relations 

<5  = | e1a " 1 1 X// " , e2  = M1/2<5y,  Xj  = 6ylt  x£  = <52|a|1/2y2, 


For  (u,<5)  in  a neighborhood  of  the  origin,  (£.,£_)  belongs  to 

1 «- 

a region  of  the  form  {(e^.e^:  | | < < (constant)  e^}. 

The  y^  are  assumed  to  lie  in  a bounded  set,  say  |y^|  < M.  A 
further  discussion  of  the  scaling  is  given  in  Section  6. 

After  scaling  (4.1.1)  becomes 


Yj  = y2  + 52g1(p,6,y) 

9 2 = s8»(e1)y1  ♦ wy2  - yj  ♦ 6>yjy,  ♦ <52s2(u , <S,y) 


(4.3.1) 


■M 
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— -1/2 

where  the  dot  means  differentiation  with  respect  to  t,  y = 3|a| 
and  g^(u,6,y)  B 0(1).  The  size  of  the  bounds  on  the  g.^  depends 
on  M,u,6.  We  write  t for  T from  now  on. 

The  cases  > 0 and  < 0 are  treated  separately. 


4.4.  The  Case  > 0 

With  > 0,  equation  (4.3.1)  becomes 


yx  = y2  + Cv» , 6 ,y) 

y2  = yx  + py2  • y\  + syyiy2  + 62g2(^»6*y) 


(4.4.1) 


Let  H(ylty2)  = (y2/2) 
of  (4.4.1), 


(y|/2)  + (y^/4) . Then  along  solutions 


H(y1,y2)  - yy2  ♦ 6y y\y\  ♦ o(62) 


(4.4.2) 


Note  that  for  u = 6 = 0,  H is  a first  integral  of  (4.4.1). 

The  level  curves  of  H(y^,y2)  = b consist  of  a figure  of  eight 

if  b * 0,  and  a single  closed  curve  if  b > 0 (see  Figure  1). 

For  b > 0,  the  curve  H(y1,y2)  = b passes  through  the  point 

1/2 

^1  = ^2  * ' * For  b > 0 and  6 sufficiently  small,  we 

? 

prove  the  existence  of  a function  u = ^(b.S)  = -yP(b)6  + 0(6  ). 

For  b > 0,  (4.4.1)  with  M = Mj(b,6)  has  a periodic  solution 

1/2 

passing  through  the  point  * 0,  y2  ■ (2b)  ' . With  u = Vi  j (0 , 6) , 
(4.4.1)  has  a figure  of  eight  solution. 

For  fixed  u,6,  the  number  of  periodic  solutions  of  (4.4.1), 
surrounding  all  three  fixed  points,  depends  upon  the  number  of 


solutions  of 


U = WjCb.6)  = -YP(b)6  + 0(6*).  (4.4.3) 

We  prove  that  P(b)  •*  <*>  as  b -*■  00  and  that  there  exists  b^  > 0 
such  that  P'(b)  < 0 for  b < b^  and  P'(b)  > 0 for  b > b^. 

These  properties  of  P(b)  determine  the  number  of  solutions  of 
(4.4.3). 

Suppose,  for  simplicity  of  exposition,  that  y^(b,6)  = -YP(b)6 
and  that  Y < 0.  If  0 < b.,  < b^,  then  there  exists  b7  > b^  such 
that  Uj(b-,,<5)  = Uj(b3>6).  Hence,  if  U = Uj(b.,,6),  then  (4.4.1) 

1/2 

has  two  periodic  solutions,  one  passing  through  = 0,  y,,  = (2b-,)  , 

1/2 

the  other  passing  through  y^  = 0,  y?  = ( 2b^ ) . If  u > 1^(0, 6), 

then  (4.4.1)  has  one  periodic  solution  surrounding  all  three  fixed 
points.  Finally,  if  u = P(bj,6),  then  the  periodic  solutions 
coincide. 

In  Figure  4,  the  periodic  solutions  surrounding  all  three  fixed 
points  in  regions  3-5  correspond  to  the  periodic  solutions  of  (4.4.1) 
which  are  parametrized  by  u = U^(b,6),  b > b^.  Similarly  the 
’’inner"  periodic  solutions  in  region  5 are  parametrized  by 
u = M1(b,6),  0 < b < bj.  The  curve  in  (e^^)  space  corre- 

sponds to  the  curve  u = p.(0,<S).  Similarly  the  curve  L, 
corresponds  to  the  curve  y = y^(b^,6)  (see  Figure  2). 

In  general  y^(b,6)  is  not  identically  equal  to  ->P(b)6, 
but  the  results  are  qualitatively  the  same.  For  example,  we  prove  the 
existence  of  a function  bj(6)  = bj  + 0(6),  such  that  if  y,6 
satisfy  y = y ^ (b^  (6)  , <5)  , then  equation  (4.4.3)  has  exactly  one 
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i . 

solution.  The  curve  y = Mj(bj(6),<5)  which  is  mapped  into  the 
curve  in  (ej,e9)  space,  corresponds  to  the  points  where  the 

two  periodic  solutions  coincide. 

- 

If  b < 0,  then  the  set  of  points  for  which  H(yj,y.,)  = b 
consists  of  two  closed  curves  surrounding  the  points  (-1,0)  and 
(1,0).  For  0 < c < 1 , we  prove  the  existence  of  a function 
6 = U 2 (c  * ^ ) = "YQ(c)6  + 0(6^).  For  y = y0(c,<$),  (4.4.1)  has  a 
periodic  solution  surrounding  the  point  (1,0)  and  passing  through 
y i = c,  y^,  = 0.  Using  f(x)  = -f(-x),  this  proves  the  existence 
of  a periodic  solution  surrounding  the  point  (-1,0)  and  passing 
through  yj  = -c,  y,  = 0. 

We  also  prove  that  Q'(c)  > 0 for  0 < c < 1.  Let  6 > 0 

and  suppose  y satisfies 

Li 

Vi  ->  C o , 6 ) < -sgn(Y)y  < U(1,<S)  . (4.4.4) 

l . L 

Then  the  equation  y = u 2 C c » <5 ) has  exactly  one  solution.  Hence, 
for  fixed  y , <5  satisfying  (4.4.4),  equation  (4.4.1)  has  exactly 
one  periodic  solution  surrounding  (1,0).  The  region  in  (y , 6) 
space,  corresponding  to  (4.4.4),  is  mapped  into  region  4 in  (l.,l.,) 
space  (see  Figure  2). 

Lemma  1.  For  6 sufficiently  small,  there  exists  a function 
y = y(<5)  = (4/5)6  ♦ 0(6")  such  that  when  y = y(6),  (4.4.1)  has 
a homoclinic  orbit. 

u 

t 

K 


Proof . Let  S(u,<5),  l)(y,6)  be  the  stable  and  unstable  manifolds  of 
the  fixed  point  (0,0)  in  (4.4.1).  These  manifolds  exist,  since 
for  m , <S  sufficiently  small,  (0,0)  is  a saddle  [231.  Let 
H(u,6,  + ) be  the  value  of  H(yj,y,)  when  U(u,<S)  hits  = 0, 

v^  > 0.  Similarly,  H(p,6,-)  is  the  value  of  H(Vj,y.,)  when 
S(y,<5)  hits  y-,  = 0,  y^  > 0.  H(M,6,±)  are  well  defined  since 

stable  and  unstable  manifolds  depend  continuously  on  parameters. 

Let  I(u,<5,  + ) denote  the  integral  of  H(y^,y.,)  over  the 
portion  of  U(u,<5)  with  v^  > 0,  y,  > 0 from  = Y2  = ^ to 
y,  = 0,  y.j  > 0.  Then 


H(m,6,-)  = I (u  , 5 , +)  . 


(4.4.5) 


Similarly,  I(u,6,-)  denotes  the  integral  of  H(yj,y,)  over  the 

portion  of  S(y,6)  with  y.  > 0,  y~  < 0 from  y.  = y,,  = 0 to 

1 Z 1 z 

>’l  = o,  y 2 > 0,  so  that  H(p,S,-)  = I(y,6,-)- 

Equation  (4.4.1)  has  a homoclinic  orbit  (with  y ^ > 0)  if  and 
only  if 


1I(U,6,  + ) - H (y  , 6 , - ) = 0. 


We  solve  (4. 4. ft)  by  the  implicit  function  theorem 
Using  (4.4.2)  and  (4.4.5), 


H(u,<5,  + ) = (uy;  + y<Sy“y“)dt  ♦ 0(y“  ♦ 6"), 

+ “ 1 “ 


(4.4.6) 


(4.4.7) 


where  the  above  integral  is  taken  over  the  portion  of  U(0,0)  from 


Yl  = >’2  = 0 to  y = 0,  y,  > 0.  Similarly, 
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L 


L. 


L 


L 

t 

l 


h 


»(»,«,-) 


[ (py2  ♦ yfiy2y2)dt  ♦ 0(p2  ♦ 62)  , 


(4.4.8) 


where  the  integral  is  taken  over  the  portion  of  S(0,0)  from 

yl  = y2  = ° t0  yi  = °»  y 2 > °*  Using  (4.4.7),  (4.4.8)  and 
U(0,0)  = -SCO, 0) , we  obtain 


H(M  , <5 , + ) = - H ( M , 6 , - ) ♦ 0 (p  2 ♦ 62) 


(4.4.9) 


Using  (4.4.7)  and  (4.4.9), 


Ip  (Hj  (0,0,-)  - Hj  (0 , 0 , -)  ) - 2 | y2dt  > 0, 


so  that  by  the  implicit  function  theorem,  we  can  solve  (4.4.6) 
for  p as  a function  of  6,  say  p = p(6).  We  now  show  how  to  get 
an  approximate  formula  for  p(<5).  We  can  write  equation  (4.4.6) 
in  the  form 


p | y2dt  + y6 


y2y^dt  ♦ 0 (p  2 ♦ 6 2 ) = 0. 


Hence,  using  (4.4.1)  and  y2  = ♦ [yj  - (yj/2)]1/2,  we  obtain 


- Y6 


>'ly2dy2 


u = 


f.y2dy 1 


♦ 0(6  2 ) = - ♦ 0(62) 


This  completes  the  proof  of  Lemma  1. 


*■ 

■■i 


1 


I 


j 
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Lemma  1 proves  the  existence  of  a homoclinic  orbit  of  (4.1.1) 
when  ei»e2  on  t^e  curve  given  by 

e2  - -(4/5)|a|'1ge1  ♦ 0(e*/2). 

Using  f(x,£)  * -f(-x,e),  when  lie  on  L^  , equation  (4.1.1) 

has  a figure  of  eight  solution. 

We  now  prove  the  existence  of  periodic  solutions  of  (4.4.1) 

surrounding  all  three  fixed  points.  In  the  introduction  to  this 

section,  we  stated  that  (4.4.1)  has  a periodic  solution  passing 

1/2 

through  = 0,  y^  = (2b)  for  any  b > 0.  In  Lemma  2,  we  only 

prove  this  for  "moderate"  values  of  b.  The  reason  for  this  is  that 
in  (4.3.1)  the  g^  are  bounded  oiuy  for  y in  a bounded  set.  In 
Section  6,  we  show  that  by  a simple  modification  of  the  scaling,  we 
can  extend  these  results  to  all  b > 0. 

Lemma  2.  Fix  > 0.  Then  for  0 < b < b and  6 sufficiently 

small,  there  exists  a function  u = Ljtb.S)  = -YP(b)6  ♦ 0(6“')  . If 
V*  = Wj(b,6)  in  (4.4.1),  then  (4.4.1)  has  a periodic  solution 
passing  through  yl  = 0,  y2  = ( 2b ) 1 ^ 2 . As  b ■*  0 the  periodic 
solution  tends  to  the  figure  of  eight  solution  obtained  in  Lemma  1. 

Proof.  Let  H(v , <5  ,b , ♦)  be  the  value  of  Hfy,  .y„)  when  the  orbit 

1/2 

of  (4.4.1),  which  starts  at  yj  = 0,  y2  = (2b)  , intersects 

y,  * 0.  Similarly,  H(u,<5,b,-)  is  the  value  of  H(y1,y,)  when  the 

1/7 

orbit  of  (4.4.1),  which  starts  at  y.  = 0,  y.,  = -(2b)  ' , is  integrated 

i.  *- 

{ 

I 
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I 


backwards  in  time  until  it  intersects  y2  = 0.  Then  (4.4.1)  has  a 
periodic  solution  passing  through  y = 0,  y2  = (2b) 1/2  if  and  only 


H(M.«,b,  + ) - H (M  , 6 , b , - ) = 0. 


(4.4.10) 


Let  l(lJ,(5»b,  + ) denote  the  integral  of  H(y^,y2)  over  the 

portion  ol  the  orbit  of  (4.4.1)  with  y.,  > 0,  starting  at 

1/2 

yl  ~ °»  y2  = (2b)  ~ and  finishing  at  v?  = 0,  y1  > 0.  Similarly, 

I (u , 6 , b , - ) is  defined  by  integrating  backwards  in  time.  Thus, 


H(P,  6,b,  ±)  = b ♦ I (u,  6, b , +)  . 


(4.4.11) 


Using  (4.4.2)  and  (4.4.11), 


,<5,b,+)  = b + (py;  + YSy2y2)dt  + 0(p2  + 62) 


(4.4.12) 


where  the  above  integral  is  taken  over  the  portion  of  the  orbit  of 
(4.4.1)  with  u = <5  = 0 from  yJ  = 0,  y2  = (2b)1//2  to  y1  = c, 
y7  = 0 where 


4b  = c4  - 2c2. 


(4.4.13) 


Similarly,  H(u,«,b,-)  = -H(u,6,b,+)  ♦ 0(m2  + <52) , so  that 
equation  (4.4.10  ) may  be  written  in  the  form 
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.2. .2 


(uy 2 + YSy^ypdt  ♦ 0(U  ♦ 6 ) = 0.  (4.4.14) 


Hence  by  the  implicit  function  theorem,  we  can  solve  (4.4.14)  to 
obtain  u = -YP(b)6  + 0(6^)  where 


P(b)  = 


I yly2dy 


y2iy  1 


(4.4.15) 


In  order  to  prove  that  the  periodic  solution  tends  to  the 
figure  of  eight  solution  as  b -*•  0,  we  prove  that 

H(w,6,b,±)  - H(p,<5,±)  as  b - 0 . (4.4.16) 

This  does  not  follow  from  continuous  dependence  of  solutions  on 
initial  conditions,  since  as  b -*■  0 the  period  of  the  periodic 
solution  tends  to  infinity.  The  same  problem  occurs  in  Kopell 

and  Howard  [39,  p.  339]  and  we  outline  their  method.  For  y.^  and 
y,  small,  solutions  of  (4.4.1)  behave  like  solutions  of  the  linear- 
ized equations.  The  proof  of  (4.4.16)  follows  from  the  fact  that  the 
periodic  solution  stays  close  to  the  solution  of  the  linearized 
equation  for  the  part  of  the  solution  with  yj»>'2  small  and  con- 
tinuous dependence  on  initial  data  for  the  rest  of  the  solution. 


Lemma  3.  P(b)  -*•  00  as  b -*■  00 . 


n 
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Proof.  The  integrals  in  (4.4.15)  are  taken  over  the  curve 
y 2 = lyj  - (y4/2)  + 2b] 1 ^ , from  y^  = 0 to  y^  = c where 
defined  by  (4.4.13).  Thus,  Jg(b)P(b)  = J^(b)  where 


'.(b)  = [\2l(w2  - (w4/2)  + 2b)1/2dw, 
1 J n 


(4.4.17) 


Substituting  w = cz  in  (4.4.17)  we  obtain 


J.  (b)  = c' 

l 


(cz)  *’1g(z)dz 


where  g(z)  = [ ( z ^ - 1 ) + (c  /2)(l-z  )] 


? - 4, ,1/2 


-1. 


Since  g(z)  < g(c  ) 


for  0 < z < 1 , we  have  that  Jg(b)  - DjC  for  some  positive  con- 


stant . Similarly,  there  exists  a positive  constant  D-,  such 


that  J^(b)  > D^c  . The  result  now  follows. 


Lemma  4 . There  exists  bj  > 0 such  that  P ’ ( b ) < 0 for  b < bj 


and  P * ( b ) > 0 for  b > b^. 


i 


1 


f « 


Proof . It  is  easv  to  show  that 


as  b -*■  0 . Hence,  by 


Lemma  3 it  is  sufficient  to  show  that  if 


0 then 


Differentiating  (4.4.17) 


with  respect  to  b we  obtain 


Integrating  by  parts  in  J.  we  obtain 


.7; 
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j . [C  iw_4.-w2)  Jw> 
0 J o rTw) 


(4.4.19) 


Also 


C liiwl  dw  . fC 
o rM  0 


* -2bl  dw.  (4.4.20) 


Similarly, 


re  ,6  .4 


3Ji  - 


d“  * 3 In  rM  l»2  - * 2b!d“ 


(4.4.21) 


Using  (4.4.18)  - (4.4.21)  we  can  express  JQ  and  J-^  in  terms 
of  Jq  and  J^.  A straightforward  calculation  yields 


3J0  = 4bJ'  ♦ J’ 


(4.4.22) 


15J1  = 4bJp  + C4+12b)J' 

Suppose  that  P • (bj_)  = 0.  Then  JQ(b1)P"(b1)  = J'1'(b1)  - 
P(b1)J,Q(b1).  Using  (4.4.22)  we  obtain 

4b1(4b1+l)[J’1’(b1)  - Ptb^J^bp]  = J^(bx)  [P2(b1)  + SbjPtbj)  - 4bj].  (4.4.23) 

Hence  P"(b..)  has  the  same  sign  as 


P (bj)  + 8b1P(b1)  - 4bj 


(4.4.24) 
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Since  P'(b^)  = 0 we  have  that  J|(b^)  = P(b^)JQ(b^).  Using 
(4.4.22)  we  obtain 

SP2(b1)  + SbjPCbj)  - 4P(b1)  - 4bj  = 0.  (4.4.25) 

Using  b^  > 0,  it  is  easy  to  show  that  (4.4.25)  implies  that 
P^)  < 1.  Using  (4.4.24)  and  (4.4.25),  P"^)  has  the  same  sign 
as  P ( b ) - P2(b1).  This  proves  that  P"^)  > 0 as  required. 

Lemma  5.  For  6 sufficiently  small  there  exist  b^(6)  = b^  + 0(6), 
b7 ( 6)  = b2  + 0(6),  where  P(0)  = P(b2),  with  the  following  properties: 
Let  6 > 0. 

(i)  If  u = u . (bj ( 6) , 6)  , then  the  equation 


U = M x (b , 6) 


(4.4.26) 


has  exactly  one  solution. 

(ii)  If  Y < 0 and  M1(b1(6),6)  < u < M1(b2(6),6),  then 
equation  (4.4.26)  has  exactly  two  solutions  b^(6),  b^(6)  with 
b^(6)  = b^  + 0(6),  b^(6)  = b^  + 0(6),  where  b^  and  b^  are 
solutions  of 


P (b)  = -u  6* 1 Y 1 


(4.4.27) 


r result  holds  if  Y 0. 
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(iii)  If  Y < 0 and  u > u ^ (b  ( 6)  , 6)  , then  (4.4.26)  has 

exactly  one  solution  b^(6)  = b^.  + 0(6),  where  b^  is  the  unique 
solution  of  (4.4.27). 

(iv)  If  Y < 0 and  u < Uj(bj(6),6),  then  (4.4.27)  has  no 

solutions.  A similar  result  holds  for  Y > 0. 

Proof . By  Lemma  4,  there  exists  b,  > 0 such  that  P(b7)  = P(0) 

and  P ' ( b T ) > 0.  Set  g(z,6)  = 6 1 (q  1 (b.,  + z , 6)  - ^(0,6)1, 

for  6*0  and  g(z,0)  = 0.  Then  g(z,6)  = P'fb^lz  + 0 1 1 6 1 + z“) . 

By  the  implicit  function  theorem  there  exists  z(6)  = 0(6)  such  that 
g(z(6),6)  = 0.  Hence,  if  b7(6)  = b , + z(<5)  then  u j ( b ( <5)  , 6)  = 

P j ( 0 , 6) . The  existence  of  b^(6)  is  proved  in  a similar  way.  The 
rest  of  the  Lemma  follows  from  the  properties  of  P(b) . 

We  now  prove  the  existence  of  periodic  solutions  of  (4.4.1) 
surrounding  a single  fixed  point. 

Lemma  6.  Por  0 < c < 1 and  6 sufficiently  small,  there  exists 
u ■ u , (c , 6)  = - YQ (c ) 6 * 0(62).  If  y = u,(c,6),  then  (4.4.1)  has  a 

periodic  solution  passing  through  y^  = c,  y,  = 0.  As  c -*■  0 the 
periodic  orbits  tend  to  the  homoclinic  orbit  obtained  in  Lemma  1. 

Proof.  Let  H(u,6,c,  + ) be  the  value  of  H(Vj,y-,)  when  the  orbit 
of  (4.4.1)  starting  at  y,  = c,  y,  « 0 intersects  y-,  c 0 , y , > 0. 
Similarly,  H(u,6,c,-)  is  the  value  of  HCy^.y-,)  when  the  orbit 
of  (4.4.1)  starting  at  yj  = c,  y,  = 0 is  integrated  backwards 
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in  time  until  it  intersects  y^  = 0.  Equation  (4.4.1)  will  have  a 
periodic  orbit  passing  through  = c>  ^2  = ^ if  and  only  if 

H (u  , 6 , c , + ) - H (p  , 6 , c , - ) = 0.  (4.4.28) 

Using  the  same  method  as  in  Lemma  2,  we  can  rewrite  equation 
(4.4.28)  as 

G (p  , 6 , c)  = | (py2  + y6y2y2)dt  + 0(p2  + 62)  = 0 (4.4.29) 

where  the  above  integral  is  taken  over  the  curve 

y2  = r(yx)  = [y2  - (yJ/2)  + (c4/2)  - c2]1/2,  (4.4.30) 

from  y^  = c,  y2  = 0 to  y?  = 0 again.  By  (4.4.29) 


fiT  G(0,0,c) 


y,dt  > 0, 


(4.4.31) 


Thus,  for  fixed  c we  can  solve  (4.4.29).  We  cannot  solve  (4.4.29) 
uniformly  in  c however,  since  as  c -*■  1 the  right  hand  side  of 
(4.4.31)  tends  to  0.  We  use  a method  similar  to  that  used  by 
Kopell  and  Howard  in  [39,  p.  337-  338]  to  obtain  ^(c.^)  for 
0 < c < 1. 

Equation  (4.4.1)  has  a fixed  point  at  (y^,y2)  = (1,0)  + 
0(|m|,|6|).  Eor  p , 6 sufficiently  small  we  make  a change  of 
variables  y^  = hj (y^ , y^ ,u , 6)  so  that  the  fixed  point  is  trans- 
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formed  into  (y1,y2)  = (1,0).  An  easy  calculation  shows  that  if 
we  make  this  transformation,  then  the  only  change  in  (4.4.1)  is 
in  the  functions  g^  and  g2-  We  suppose  that  the  above  change 
of  variables  has  been  made  and  we  write  >'j  for  y^. 

The  curves  H(y1,y2)  = H(c,0)  can  be  written  in  the  form 

4H(y1,y2)  ♦ 1 - 2y2  + (yx  - 1) 2 (y^l)  2 = (c-l)2(c+l)2.  (4.4.32) 

Thus,  for  c close  to  1 the  closed  curves  are  approximately 


and  - 1.  Hence  there  exist  functions  q-(M,«S)  = 0(|p|  + | <5 1 ) 

such  that 

qi(M,<5)  - 3T  ln(4HCy1,y2)  + i)  5 q20,6). 

Integrating  the  above  inequality  over  the  curve  given  by  (4.4.30) 
we  obtain 


exp(qi(p,6)T)  < (4H(M,c,+)  + l]/[4H(c,0)  + 1]  < exp(q2(p , 6)T)  (4.4.35) 


where  T is  a bound  for  the  time  taken  to  trace  the  orbit.  Using 
(4.4.29),  (4.4.32)  and  (4.3.35)  we  see  that  G(u,6,c)  is  bounded. 
The  fact  that  ( 9/ 9u ) G (c , 0 , 0)  is  bounded  away  from  zero  follows 
easily  from  (4.4.32).  This  proves  that  for  6 sufficiently  small 
and  0 < c < 1 , we  can  solve  (4.4.32)  to  obtain  y = m2(c,6)  = 
-YQ(c)6  + 0(62)  where  Jq(c)Q(c)  = Jjfc), 


JjCc) 


d 

[ w^r(w)dw 


and  r(w)  is  defined  by  (4.4.30)  and  r(d)  = 0,  d > c. 

The  fact  that  the  periodic  solution  tends  to  the  homoclinic 
orbit  is  proved  in  the  same  way  as  the  corresponding  result  in  Lemma 
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Using  f (x)  = -f(-x).  Lemma  6 proves  the  existence  of  periodic 
solutions  of  (4.4.1)  surrounding  (-1,0). 

Lemma  7 . Q'(c)  > 0 for  0 < c < 1. 

Proof . We  write  Q»Jq  and  as  functions  of  b where 

4 2 

4b  = c - 2c  . Since  (db/dc)  < 0 for  0 < c < 1 , we  must  prove 
that  Q'(b)  < 0 for  -1  < 4b  < 0.  Following  the  same  procedure 
as  in  Lemma  4,  we  find  that  JQ,J^,  Jl,Jl  satisfy  equation  (4.4.22). 
Thus,  if  Q^bj)  = 0 then 

SQ2^)  + 8b1Q(b1)  - 4Q(b1)  - 4bx  = 0.  (4.4.36) 

Since  -1  < 4bj  < 0,  the  roots  of  (4.4.36)  are  less  than  1.  Hence, 
if  Q’^)  = 0 then  Q(bj)  < 1.  Also,  from  (4.4.23),  if  Q'Cbj)  = 0 
then  Q"(bj)  has  the  same  sign  as  4bj  - 8b1Q(b1)  - Q2(bj).  Using 
(4.4.36)  and  Q ( b ^ ) < 1,  this  implies  that  Q’^bj)  < 0.  Since 
Q ( - 1 / 4 ) = 1,  Q ( 0 ) = 4/5,  this  shows  that  Q’(b)  < 0 for  -1  < 4b  < 0. 
This  completes  the  proof  of  the  Lemma. 

4.5.  The  Case  < 0 

With  < 0,  equation  (4.3.1)  becomes 

y1  ■ y2  * 

y 2 = *yi  + py2  ■ yj  + 6yy2y2  + 62s2k  »5»y) . 


(4.5.1) 
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Let  H1Cy1,y2)  = (y2/2)  + (y^/2)  + (y*/4) . Then  along  solutions 
of  (4.5.1),  Hj  = uy2  + ^Yy^y^  + Using  the  same  methods  as 

in  the  previous  Lemmas,  we  prove  that  (4.5.1)  has  a periodic  solu- 
tion passing  through  y^^  = c > 0,  y2  = 0,  if  and  only  if 
U = U j (c , 6)  = - yR(c)  6 + 0(62),  where  Jfl(c)R(c)  = ^(c), 


Jj(c) 


fc  2 i 

w^  r (w)dw, 

J0 


r(w)  = [2b  - (w^/2)  - w2]*^2,  4b  = c^  + 2c2. 


In  order  to  prove  that  for  fixed  p,6,  equation  (4.5.1)  has 
at  most  one  periodic  solution  we  prove  that  R is  strictly 
monotonic . 


Lemma  8 . R' (c)  > 0 for  c > 0. 

Proof . We  write  R,Jq,J^  as  functions  of  b.  It  is  sufficient 
to  prove  that  R'(b)  > 0 for  b > 0. 

Using  the  same  methods  as  before,  we  show  that 


3Jq  = 4bJ*  - JJ 


lSJj  = (12b  + 4)J’  - 4bJ p 


(4.5.2) 


3v/dw  = 4bJ,  . 4J, 
r (w)  0 1 


0 


(4.5.3) 


Now  R'  has  the  same  sign  as  S where  S = 15  [Jj . By 
(4.5.2), 

S = (8b-4) JqJ|  = 5(J{)2  + 4b  (Jq)  2 . 

By  (4.5.3),  bJg  > JJ.  Hence,  for  b > 2. 

S > (JJ)2b_1(8b-4-5b+4)  = 3(Jp2  > 0. 

2 2 

Similarly,  if  0 < b < 2,  then  S > 3b  (Jq)  > 0.  This  completes 
the  proof  of  Lemma  8. 

4.6.  More  Scaling 

In  Section  4 we  proved  the  existence  of  periodic  solutions 

1/2 

of  (4.4.1)  which  pass  through  the  point  (0,(2b)  ).  Lemma  3 

indicates  that  we  may  take  b to  be  as  large  as  we  please.  How- 
ever our  analysis  relied  on  the  fact  that  the  g^  are  0(1) 
and  this  is  true  only  for  y in  a bounded  set.  Also,  our  analysis 
in  Sections  4 and  5 restricts  ei>e2  t0  a reRi°n  °f  the  form 
{(el»e2):  I e1 I 5 e0>  e1  5 (constant)  e^}. 

To  remedy  this  we  modify  the  scaling  by  setting 
6 = |e,<*  *|^2h  where  h is  a new  parameter  with  0 < 4h  < 1 
say.  The  other  changes  of  variables  remain  the  same.  Note  that  if 
U,6  lie  in  a full  neighborhood  of  the  origin  then  (ei*e2^  anc* 
(x^,X2)  lie  in  a full  neighborhood  of  the  origin. 

After  scaling,  (4.1.1)  becomes 
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Y1  = y2  + <s2g1Oi»5,h,y) 

y2  = h2sgn(e1)y1  + Uy2  - 


(4.6.1) 

yl  + 6Yyly2  + 62g2(M> 6»h*y) • 


For  0 < 4h  < 1 and  y,6  sufficiently  small,  the  g^  are  0(1) 
for  y in  a bounded  set. 

Let  ex  > o.  Let  H^.y^  = (y2/2)  - (h2y2/2)  + (yJ/4)  . 

Then  along  solutions  of  (4.6.1),  H2  = M y2  + y<sy2Y2  + 0(M2+<52). 
Following  the  same  procedure  as  in  Section  4,  we  find  that  (4.6.1) 
has  a periodic  solution  passing  through  y1  = 0,  y2  = 1,  if  and 
only  if 

w = h2M1(b,6)  = - yh2P (b) 6 + 0(62), 

where  2b  = h 4 . An  easy  computation  shows  that  as  h -*•  0, 
h2P(b)  - K ♦ 0 (h2 ) , where  JQK  = /I 

J.  = f 1w2l(l-w4)1/2dw. 

1 J 0 

Thus,  for  small  h,  (4.6.1)  has  a periodic  solution  passing  through 
yl  = 0,  y2  = 1 if  and  only  if  y = -yK6  + 0 (h 2 1 6 1 + 62) . In 
particular,  when  e = °,  e2  > °,  & < 0,  (4.1.1)  has  a periodic 
solution  passing  through  x 1 = 0,  x2  = |u|(-gK)’1e2  + Ofe^2). 

For  < 0,  a similar  anlysis  shows  that  (4.6.1)  has  a periodic 
solution  passing  through  yx  = 1,  y2  = 0,  if  and  only  if 

y = h2y  2 (h* 1 , 6)  = -yh2Q(h'1)6  + 0(62), 


[ 
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c'-T~ 


and  that  as  h -*•  0,  h“Q(h  *)  = (K//2)  + 0(h"). 

4.7.  Completion  of  the  Phase  Portraits 

Our  analysis  in  the  previous  sections  proves  the  existence  of 
periodic  and  homoclinic  orbits  of  (4.1.1)  for  certain  regions  in 
£1  ' *"2  sPace-  We  now  show  how  to  complete  the  pictures. 

Obtaining  the  complete  phase  portrait  in  the  different  regions  in- 
volves many  calculations.  However,  the  method  is  the  same  in  each 
case  so  wc  only  give  one  representative  example.  We  prove  that  if 
Y < 0,  then  the  "outer"  periodic  orbits  in  region  5 is  stable.  We 
use  the  scaling  given  in  Section  3. 

Since  we  are  in  region  5,  ^(0,6)  > u > y ^ (bj  ( 6)  , <5) . Fix  y 
and  6 and  let  b-,  be  the  solution  of  y = y^(b,6)  with 

b-,  > b ^ ( 6 ) . Then  we  have  to  prove  that  the  periodic  solution  T 

1/  ? 

of  (4.4.1)  passing  through  y ^ = 0,  y,  = (2b-,)  ' is  stable. 

Let  b-  be  the  solution  of  y = y^(b,6)  with  b,  < b^(6). 

Then  there  is  a periodic  solution  passing  through  Vj  = 0, 

1/2 

y-,  = ( 2b ^ ) . In  fact,  we  prove  that  any  solution  of  (4.4.1) 

starting  "outside"  this  periodic  solution  (and  inside  some  bounded 
set)  tends  to  r as  t . 

Let  b > bj.  Let  c(b,+)  be  the  value  of  y.  when  the  orbit 

of  (4.4.1)  starting  at  y^  * 0,  y,  = (2b)*^“  first  hits  y-,  = 0. 

Similarly,  c(b,-)  is  the  value  of  y ^ when  the  orbit  of  (4.1.1) 

1/2 

starting  at  y^  = 0,  y,  * (2b)  is  integrated  backwards  until  it 

hits  y 2 = 0.  Using  f(x,e)  = -f(-x,e),  the  solution  passing  through 
(0,(2b)A/4),  spirals  inwards  or  outwards  according  to  the  sign  of 


I 

I 

I 

! 


J 


i 1 


I 
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c(b,+)  + c(b,-).  Using  the  calculation  in  Lemma  2,  c(b,+)  + cCb,-) 
has  the  same  sign  as  H(p  , <5 , b , +)  - H(p,<5,b,-)  which  in  turn  has 
the  same  sign  as 


S = (b2, 6)  - M x (b , 6) . 

Using  the  properties  of  u^(b,6)  given  in  Lemma  5,  S is  positive 
if  b < b^  and  negative  if  b > b2-  The  result  now  follows. 

4.8.  Remarks  and  Exercises 

Remark  1.  Consider  equation  (4.1.1)  under  the  hypothesis  that  the 
linearization  is  given  by  x = A(e)x  where  A(e)  is  defined  by 
(4.1.2).  Suppose  also  that  (4.2.1)  and  (HI)  hold.  Make  the  change 
of  variables  x = B(c)x,  where  B(e)  = I - rct^Afe)  and 

a3f1(0,0)  93f?(0,0j 

r = 3 , a = ^ . 

3xJ  3xJ 

The  map  x ■+  x will  be  one-to-one  for  c sufficiently  small. 
Using  the  fact  that  A(c)  and  B(e)  commute,  it  is  easy  to  show 
that  the  transformed  equation  satisfies  all  the  above  hypotheses 
and  that  in  addition  it  enjoys  ( H 3 ) . In  particular,  if  the  trans- 
formed equation  is 


x = F(x,e) 


then 
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33F1 CO, 0) 


93F2(0,0)  3 3f 2 (0 , 0) 

3x3  3x3 

33F2(0,0)  33f2(0,0)  333f1(0,0) 

3^3x2  ~Vx“3x2  3x^ 

Reraark_2.  We  have  assumed  that  f(x,e)  = -f(-x,e).  If  we  assume 
that  this  is  true  only  for  the  low  order  terms  then  we  would  obtain 
similar  results.  For  example,  on  Lx  we  would  get  a homoclinic 
orbit  with  x^  > 0.  Similarly  on  another  curve  LJ  we  would  get 
a homocl mic  orbit  with  Xj  < 0.  In  general  and  Lj  would  be 

different  although  they  would  have  the  same  linear  approximation 
U = (4/5)6. 

Remark_3.  Suppose  that  we  only  assume  (HI)  and  (H3) . Then  we  can 
still  obtain  partial  results  about  the  local  behavior  of  solutions. 
For  example,  in  Lemma  2 we  did  not  use  the  hypothesis  3 nonzero. 
Hence,  for  each  b > 0,  (4.1.1)  has  a periodic  solution  through 
X1  = °»  x2  = M 1//2(2b)1/2ex  for  some  Ej  and  e2  with  Ej  > 0. 
However,  we  cannot  say  anything  about  the  stability  of  the  periodic 
orbits  and  we  cannot  say  how  many  periodic  orbits  (4.1.1)  has  for 
fixed  ex  and  e2« 
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Exercises . 

(1)  Suppose  that  and  8 are  negative.  Prove  that  for 
(£1,£2^  in  re8i°ns  5 and  6»  (4.1.1)  has  a connecting  orbit. 

Hint : Use  the  calculations  in  Lemma  1 and  Lemma  6. 

(2)  Suppose  that  a and  6 are  negative  and  that  (Ej,^) 
is  in  region  3.  Let  U be  the  unstable  manifold  of  the  point 
(0,0)  in  (4.1.1).  Prove  that  U is  in  the  region  of  attraction 
of  the  periodic  orbit. 

(3)  Suppose  that  u is  positive  and  8 is  negative.  Show 
that  the  bifurcation  set  and  the  corresponding  phase  portraits  are 
as  given  in  Figures  3 and  5. 


4.9.  Quadratic  Nonlinearities 

In  this  section  we  discuss  the  local  behavior  of  solutions  of 
(4.1.1)  when  the  nonlinearities  are  quadratic.  Most  of  the  material 
in  this  section  can  be  found  in  [39] . 

Suppose  that  the  linearization  of  (4.1.1)  is  x = A(e)x  where 
A(t)  is  given  by  (4.1.2)  and  that 


f(0,e)  = 0, 


fj(0,0)  = 0. 


1-  (0,0)  * 0, 


We  also  assume  that  g > 0,  u < 0;  the  results  for  the  other  cases 
are  obtained  by  making  use  of  the  change  of  variables,  t •*  -t, 


-*  -e £ » xi  -*•  ±X|  , x?  -*•  +x 


1*  2 
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Introduce  parameters  p,6,  scaled  variables  y^,y2  and  a new 
time  t by  the  relations 


6 


Cla 


-1,1/2 


x1  = 6“y1,  x2  = |a|1/“63y2> 
|a|1/26u>  t - t|a|-l/26-l. 


We  asume  that  > 0 in  what  follows,  see  Exercise  8 for  the 

case  < 0. 

After  scaling,  (4.1.1)  becomes 


yl  = y2  + 0(6“) 

y2  = Yi  + ay2  - y^  + 6Yy1y2  + 0(6“) 


(4.9.1) 


where  the  dot  means  differentiation  with  respect  to  t,  the 

1/2 

lie  in  a bounded  set  and  y = B|a|  . Note  that  by  making  a 
change  of  variable  we  can  assume  that  (1,0)  is  a fixed  point  of 
(4.9.1)  for  p,6  sufficiently  small. 

The  object  of  the  following  exercises  is  to  show  that  the 
bifurcation  set  is  given  by  Figure  6 and  that  the  associated  phase 
portraits  are  given  by  Figure  7. 

Exercises . 

(4)  Let  HCylty2)  ■ (y|/2)  - (yJ/2)  ♦ (yJ/3).  Show  that 
along  solutions  of  (4.9.1), 


H = uy“  + SYy-^y;  + 0(6“)  . 


>8  ro 


■ 
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(5)  Prove  that  there  is  a function  p = P^(c,6),  0 < c 5 1, 
such  that  if  p = p^Cc^)  in  £4.9.1)  then  there  is  a periodic 
solution  through  (c,0)  if  0 < c < 1 and  a homoclinic  orbit  if 
c = 0 . 

(6)  Prove  that  p^Cc,^)  = -yP(c)6  + 0(62)  where 
J0(c)P(c)  = J1£c), 


Ji(c)  = 


rci 


w R(w)dw, 


R(w)  = [w2  - (2/3)w3  + 2b]1/2, 

6b  = c2 (2c- 3)  , R(Cj)  = 0,  c2  > c. 


(7)  Prove  that  P(0)  = (6/7),  P (1)  = 1 and  P'(c)  > 0 for 
0 < c < 1.  Deduce  that  for  fixed  P,6,  (4.9.1)  has  at  most  one 
periodic  solution.  (To  prove  that  P'(c)  > 0 we  can  use  the  same 
techniques  as  in  the  proof  of  Lemma  4.  An  alternative  method  of  proving 
P'(c)  >0  is  given  in  (39].) 

(8)  If  < 0,  then  after  scaling  (4.1.1)  becomes 

t1  = y2  + 0 ( 62 ) 

y 2 = -yx  + py2  * y\  + vfiy^  + °C62)- 


Put  y1=z-l,y=p-y6.  Then 

z - y2  + 0(62) 

y2=  z + py2  - z2  + 6yzy2  + 0(62) 


which  has  the  same  form  as  equation  (.4.9.1)  and  hence  transform  the 
results  for  > 0 into  results  for  the  case  < 0. 

(9)  Show  that  the  bifurcation  set  and  the  corresponding  phase 
portraits  are  as  given  in  Figures  b-7. 


CHAPTER  5 

APPLICATION  TO  A PANEL  FLUTTER  PROBLEM 


5.1.  Int  roduction 

In  this  chapter  we  apply  the  results  of  Chapter  4 to  a 
particular  two  parameter  problem.  The  equations  are 


* * Ax  ♦ f(x) 


(5.1.1) 


where 


^l’VW  ’ f(x)  * If1(x).f2(x),f3(x),f4(x)lT, 


0 

al 

0 

■c 


1 

bl 

0 

0 


0 

c 

0 

a. 


0 

0 

1 

b- 


fi(x)  - 

f3(x)  ■ 0 

9 

f2(x)  . 

Xjg(x)  , 

f4(x)  = 

4x3g(x) , 

2g(x)  - 

-AkXj  ♦ 

°*lx2 

8P 

_2 . 2 

c * 

X • aj  ■ 

* -n  y 

bj - 

* l°TI4  j 4 + 

/p  6]; 

o > 0 

are  fixed 

and  p 

0.1, 


I 
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results  from  a two  mode  approximation  to  a certain  partial  differ- 
ential equation  which  describes  the  motion  of  a thin  panel. 

Holmes  and  Marsden  [32,34]  have  studied  the  above  equation  and 
first  we  briefly  describe  their  work.  By  numerical  calculations, 
they  find  that  for  P - po  “ 1°8,  r ’ ro  “ -2.23n2,  the  matrix  A 
has  two  zero  eigenvalues  and  two  eigenvalues  with  negative  real 
parts.  Then  for  |p-pq|  and  | r - r Q I small,  by  centre  manifold 
theory,  the  local  behavior  of  solutions  of  (5.1.1)  is  determined  by 
a second  order  equation  depending  on  two  parameters.  They  then 
use  some  results  of  Takens  [49]  on  generic  models  to  conjecture 
that  the  local  behavior  of  solutions  of  (5.1.1)  for  |p-Pq|  and 
| r-r0|  small  can  be  modelled  by  the  equation 

••  • . 2 • 3 

u ♦ au  ♦ bu  ♦ u u ♦ u “0 

for  a and  b small. 

Recently,  this  conjecture  has  been  proved  by  Holmes  [33],  in 
the  case  o = 0 , by  reducing  the  equation  on  the  centre  manifold 
to  Takens  normal  form.  We  use  centre  manifold  theory  and  the 
results  of  Chapter  4 to  obtain  a similar  result. 

5.2.  Reduction  to  a Second  Order  equation 

The  eigenvalues  of  A are  the  roots  of  the  equation 

X4  ♦ djX3  ♦ d2X2  ♦ d3\  ♦ d4  - 0 (5.2.1) 


j 

i 


\ 

i 


where  the  d.  are  functions  of  r and  p. 
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eiucnvalues  then  d.  = d,  = 0.  A calculation  shows  that  if  d,  = d.  = 0, 
' .->4  _•>  4 ’ 


aja2  + c = 0 


alb2  + bla2  = 0 


(5.2.2) 


or  in  terms  of  r and  P, 


47T^(7i2+r)  (4*2+r)  + p2  = o 


(5.2.3) 


(16aTi4+6p1/2)  (ti2  + D + 4(an2  + 6p1/2)  (4n2  + r)  = 0.  (5.2.4) 


We  prove  that  (5.2.3),  (5.2.4)  has  a solution  T = Tg,  p = Pg. 

From  (5.2.3)  we  can  express  P in  terms  of  r.  Substituting 
this  relation  into  (5.2.4)  we  obtain  an  equation  H(T)  * 0. 
Calculations  show  that  HfTj)  < 0,  H(T2)  > 0 where  = -(2.225)* 
and  r2  = -(2.23)*2,  so  that  H(Tg)  = 0 for  some  Tg  € (T 2 , T ^ ) . 
Further  calculations  show  that  (5.2.3),  (5.2.4)  has  a solution 
rQ,Pg  with  107.7  < Pg  < 107.8. 

In  the  subsequent  analysis,  we  have  to  determine  the  sign  of 
various  functions  of  rQ  and  Pg.  Since  we  do  not  know  Tg  and 
Pg  exactly  we  have  to  determine  the  sign  of  these  functions  for 
Tg  and  Pg  in  the  above  numerical  ranges. 

When  T = Tg,  P * Pg,  the  remaining  eigenvalues  of  A are 
given  by 


then 


( 


Li 

y 


(bj*b2)  ± [(brb2)2  ♦ 4(a^a2)]1/2 


and  a calculation  shows  that  they  have  negative  real  parts  and 
non-zero  imaginary  parts. 

We  now  find  a basis  for  the  appropriate  eigenspaces  when 
T = rQ,  P = Pq.  Solving  A Vj  = 0 we  find  that 

Vj  - ( 1 , 0 , - aj/c  , 0] T . (5.2.5) 

The  null  space  of  A is  in  fact  one-dimensional  so  the  canonical 

2 

form  of  A must  contain  a Jordan  block.  Solving  A v,  = 0 we 
obtain 

v,  = (0,1 , -bj/c  ,-aj/cl F,  AVt  = Vj.  (5.2.t>) 

The  vectors  v^  and  v,  form  a basis  for  the  generalized  eigenspace 
of  A corresponding  to  the  zero  eigenvalues.  Similarly,  we  find  a 
(real)  basis  for  the  space  V spanned  by  the  eigenvectors  correspond- 
ing to  A,  and  A ^ . Solving  Az  = A^z,  we  find  that  V is  spanned 
by  v-  and  where 

2v3  = z ♦ z,  2v4  = i(z-z)  , 

z = Il*Aj,w,AjWjr  (5.2.7) 

wc  = b,A3  - bjb2  ♦ a2. 

Let  Aq  denote  the  matrix  A when  P = P(1  and  p = p^. 
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Let  S = [v'1,v2»v3» v4  1 where  the  are  defined  by  (5.2.5), 

(5.2.6)  and  (5.2.7)  and  set  y = S *x.  Then  (5.1.1)  can  be  written 
in  the  form 


X - By  ♦ S* 1 (A- AQ) Sy  * F(y,T,P)  (5.2.8) 

where  F(y,T,P)  = S*1f(Sy), 


and  where  ^3  = Pj  + ip2»  Pi  < 0,  p2  ♦ 0. 

Then  for  | r - T q | and  j p - p Q | sufficiently  small  (5.2.8) 

has  a centre  manifold  y3  = hj (y2 ,y 2 ,T ,P) , y4  * h2 (yx ,y2 ,T , P) . The 
flow  on  the  centre  manifold  is  governed  by  an  equation  of  the  form 


yl 

0 1 

yl 

yi 

ss 

* H(r,p) 

y2 

0 

1 

y2 

y2 

+ N(y^,y2>r>P) 


(5.2.9) 


where  E(T,p)  is  a 2*2  matrix  with  E(rQ,P0)  = 0 and 
N(y1,y2,r,P)  contains  no  linear  terms  in  y^^  or  y2>  We  show 
that  there  is  a nonsingular  change  of  variables  (T,P)  (e1,e2) 

for  (T,P)  close  to  ( r q , P q ) and  a T,P  dependent  change  of 
variables  (yj,y2)  -*•  (yj^)  such  that  the  linearized  equation 
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corresponding  to  (5.2.9)  is 


(5.2.10) 


The  transformation  (Vj.y.,)  -*  Cy j ,v-,)  is  of  the  form 

Identity  ♦ 0 (.  | r - r ^ J , |P-Pq|)  and  Y = r^,  P=  p^  is  mapped  into 

fc1  = 0.  L->  = 0.  After  these  transformations,  (5.2.9)  takes  the  form 


C1  E2 


* N(y1.y2.ro»po)  * N(yi»y2’r'p)  (5.2.11) 


where  we  have  dropped  the  bars  on  the  y . . The  function  N will 
contain  no  linear  terms  in  y1  or  y,  and  N(yx , >% , rQ ,PQ)  = 0. 
Since  the  nonlinearities  in  (5.1.1)  are  cubic,  the  same  will  be 
true  of  N and  N. 

Let 


N(y1,y:,r0,P0)  = [N1(y1,y2,r0,P0),  N2(yi,y2,r0,P0)]  (5.2.12) 


and  let 


01 . 


ayl 


N2(o,o,r0,P0) 


r = 


3y 


3 


3* 


~2 — N2(o,o,r0,P0) 


s y J 3 y -> 


p = 3r  ♦ p. 


(5.2.13) 


- 


-100- 


Using  the  results  in  Chapter  4 (see,  in  particular.  Remark  1 in 
Section  8),  if  and  8 are  non-zero,  we  can  determine  the 

local  behavior  of  solutions  of  (5.2.11).  By  Theorem  2 of  Chapter  2, 
this  determines  the  local  behavior  of  solutions  of  (5.2.8). 
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of  the  Jacobian  of  the  mapping  given  by  (5.3.1).  Hence  J is 
non-zero,  so  by  the  implicit  function  theorem  we  can  use 
*-j  = -det  (C(T  , P)  ) , = t race  (.C  (T  , P) ) as  our  bifurcation  parameters. 

Approximate  formulae  for  t-j  and  t,  can  easily  be  found  if  so 
desired . 

Let  C(r,P)  = [ci ^ ) , 


M = 


0 


C11  C1 2 


L. 

-J 

r 1 

r -m 

n 

= M 

*1 

y2 

Then  the  linearized  equation  corresponding  to  (5.2.9)  is 


Note  that  M is  equal  to  the  identity  matrix  when  e.  = e7  = 0. 

1 4. 


5.4.  Calculation  of  the  Nonlinear  Terms 

We  now  calculate  the  nonlinear  term  in  (5.2.9)  when  T = T 
P = Pq.  Since  the  nonlinearities  in  (5.1.1)  are  cubic,  the  centre 
manifold  has  a "cubic  zero"  at  the  origin.  Using  x = Sv,  on  the 
centre  manifold 

—— — — — — 
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X1  ~ yl  + °3»  x2  = y2  + °1 
•a,  b. 


(5.4.1) 


1 D1  'al 

' yl  ‘ T y2  + °3»  X4  = T"  y2  + 0 


where  03  = OdyJ3  + |y2|3) 


Let  S = ttjj]  and  let 


Fj  ^yl,y2)  ~ ^yl,y2,hl^yl,y2’ r0*P0^ 


Then  using  the  notation  introduced  in  (5.2.12) 


Nity1,y2.r0-po>  * 

N2(>Wr0-p0)  = t22F2t>V>'2>  * t24F4f>'l'>’2) 


(5.4.2) 


Using  (5.4.1) 


F2(yl’y2}  = V fkyi  + 4k(-^)  2y3  + ay2y2  + 


4aa^ 


. tt^  2 ^ j 2 3 2 81c3.^b^  2 -r  _ 2 

prVz  * -rT  >V2> 


+ terms  in  + C>5 


-4a 


"T“i  2b,  n - a.  “ 

F4(yl,y2)  = ~~E~  F2(yl’y2)  + c — [kyiy2  + 4k(“T^  yiy2] 

+ terms  in  y^2  and  y3  + 05 


where  Oj.  = 0(|y^|^+  I y 2 ! ^ ) • Note  also  that  since  r 
P = PQ,  (5.2.2)  holds. 

From  (5.2.13)  and  (5.4.2) 


Tq  and 


grim'd  i •* 


1. 

L 

[ 


l 0 .N  - 


- n 


mi  * u-J-ru t22  ■ a t24) 


- 11  ^ 1 ? *•  1 

r - — kCl  ♦ 4 ( —)  j 1 1 j 2 - 4-~-  tly() 


14 


- n 


Ska  jb  j 4oa“ 


^n4 


♦ -V-)UJ2  - 4 t24)  ♦ b1kt24CcS4a“) 


Routine  calculations  show  that 


12 

as 

m“b  , ( 2a  j - b j h ,) 

1 2 

a 

ma  . 

fa 

2 

14 

a 

m c(bj»h2) 

24 

a 

-me 

m 

a 

(al+a2'bjb2)  1 

Using  (5.2.2)  and  numerical  calculations,  we  find  that 

l*  1 

t22  - 4 — t,4  = m(4rtj+a2)  > 0 


12 


al 


t,,  ' 4 c t14  " b 1 m“ C 2a  , - - 4a 1 ) > 0 


so  that  Oj  and  r are  negative.  Similarly,  the  coefficient 
of  k in  ^ is 


^2" 

c 


11 

N—  a.  b. m(  1 2a,  *a  ,)  < 0, 
~ II  I *. 


Hence  and  P arc  negative  and  the  local  behavior  of  solutions 


of  (5.1.1)  can  he  determined  using  the  results  of  Chapter  4 


< 


CHAPTER  6 

INFINITE  DIMENSIONAL  PROBLEMS 


6.1.  Introduction 

In  this  chapter  we  extend  centre  manifold  theory  to  a class 
of  infinite  dimensional  problems.  For  simplicity  we  only  consider 
equations  of  the  form 


w = Cw  + N(w) , w(0)  C Z, 

where  Z is  a Banach  space,  C is  the  generator  of  a strongly 
continuous  semigroup  on  Z and  N:  Z -*•  Z is  smooth.  (For  a 
discussion  of  semigroup  theory  see  (4,  44].)  With  appropriate 
spectral  assumptions  on  C,  it  is  very  easy  to  prove  the  analogues 
of  Theorems  1-3  of  Chapter  2.  For  generalizations  to  other 
evolution  systems  see  [30,42]. 

6.2.  Centre  Manifolds 

Let  Z be  a Banach  space  with  norm  | | *|  | • We  consider 
ordinary  differential  equations  of  the  form 

w = Cw  + N(w) , w(0)  € Z,  (6.2.1) 

where  C is  the  generator  of  a strongly  continuous  linear  semigroup 

S(t)  and  N:  Z -*■  Z has  a uniformly  continuous  second  derivative 

with  N(0)  = 0,  N'(0)  = 0 [N'  is  the  Frechet  derivative  of  N] . 

It  can  be  shown  [4,7]  that,  with  an  appropriate  definition  of  weak 


1 
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solution,  a function  w £ C([0,T];Z),  T > 0,  is  a weak  solution 
of  (6.2.1)  if  and  only  if  w satisfies 


■ t 

w(t)  = S(t)w(0)  + S(t-x)N(w(x))dt 

0 


(6.2.2) 


on  [0,T].  A standard  argument  shows  that  there  is  a unique  solu- 
tion of  (6.2.2)  defined  on  some  maximal  interval  [0,T  ),  T >0, 

m m 

and  that  if  T < <®  then 
m 


Lim  | | w(t)  | | = 00  . 

t-T" 

m 


Furthermore,  if  w(0)  £ D(A)  then  w(t)  is  differentiable  on 
(0,Tm)  and  satisfies  (6.2.1). 

To  put  (6.2.1)  into  canonical  form  we  make  some  spectral 
assumptions  about  C.  We  assume  from  now  on  that: 

(i)  Z = X ® Y where  X is  finite  dimensional  and  Y is 
closed . 

(ii)  X is  C-invariant  and  that  if  A is  the  restriction  of 
C to  X,  then  the  real  parts  of  the  eigenvalues  of  A 
are  all  zero. 

(iii)  If  U(t)  is  the  restriction  of  S(t)  to  Y,  then  Y 

is  U(t) -invariant  and  for  some  positive  constants  a,b. 


|U(t)| | < ae"D  , t > 0, 


(6.2.3) 


Let  P be  the  projection  on  X along  Y.  Let  B = (I-P)C 
and  for  x £ X,  y € Y,  let  f(x,y)  = PN(x+y),  g(x,y)  = (I-P)N(x+y). 
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liquation  (.0.2.1)  can  be  written 


x = Ax  f(x,y) 

} = By  + g(x,y)  . 


(6.2.4) 


An  invariant  manifold  for  (6.2.4)  which  is  tangent  to  X 
space  at  the  origin  is  called  a centre  manifold. 


Theorem  1 . There  exists  a centre  manifold  for  (6.2.4),  y = h(x), 

2 

| x | <6  where  h is  C . 


The  proof  of  Theorem  1 is  exactly  the  same  as  the  proof  given 
in  Chapter  2 for  the  corresponding  finite  dimensional  problem. 

The  equation  on  the  centre  manifold  is  given  by 

u = Au  ♦ f (u,h (u) ) . (6.2.5) 

In  general  if  y(0)  is  not  in  the  domain  of  B then  y(t)  will 
not  be  differentiable.  However,  on  the  centre  manifold  y(t)  = 
h(x(t)),  and  since  X is  finite  dimensional  x(t),  and  consequently 
y(t),  are  differentiable. 


Theorem  2 (a)  Suppose  that  the  zero  solution  of  (6.2.5)  is  stable 
(asymptotically  stable)  (unstable).  Then  the  zero  solution  of 
(6.2.4)  is  stable  (asymptotically  stable)  (unstable). 

(b)  Suppose  that  the  zero  solution  of  (6.2.5)  is  stable 
Let  (x(t),y(t))  be  a solution  of  (6.2.4)  with  | | (x(0) ,y (0) | | 


I 
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sufficiently  small.  Then  there  exists  a solution  u(t)  of  (6.2.5) 


such  that  as  t 


x ( t ) = u ( t ) + 0(e'Yt) 


y (t ) = h(u(t) ) + 0(e  rt) 


(6.2.6) 


where  Y > 0. 


The  proof  of  the  above  theorem  is  exactly  the  same  as  the  proof 
given  for  the  corresponding  finite  dimensional  result. 

Using  the  invariance  of  h and  proceeding  formally  we  have 

that 


h ' (x) (Ax  + f (x ,h (x) ) ] = Bh (x)  + G(x,h(x)).  (6.2.7) 


To  prove  that  equation  (6.2.7)  holds  we  must  show  that  h(x)  is  in  the 
domain  of  B. 

Let  xQ  € X be  small.  To  prove  that  h(xQ)  is  in  the  domain 
of  B it  is  sufficient  to  prove  that 


Lim 
t-0  + 


U(t)h(xn)  - h(xn) 


exists.  Let  x(t),  y(t)  = h(x(t))  be  the  solution  of  (6.2.4)  with 

x ( 0 ) = xQ.  As  we  remarked  earlier,  y(t)  is  differentiable.  From 
(6.2.4) 

y(t)  - U(t)h(x0)  + | U(t-T)g(x(T) ,y(T))dT, 
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so  it  is  sufficient  to  prove  that 

Lim  \ f u(t-T)g(x(T) ,yCT))dT 

t-0+  1 Jo 

exists.  This  easily  follows  from  the  fact  that  U(t)  is  a strongly 
continuous  semigroup  and  g is  smooth.  Hence  h(xQ)  is  in  the 
domain  of  B. 


Theorem  3.  Let  <J>  be  a map  from  a neighborhood  of  the  origin 

in  X into  Y such  that  <t>(0)  = 0,  = 0 and  <J>  (x)  € D(B) . 

Suppose  that  as  x -*•  0,  (M4>)(x)  = 0(|x|q),  q > 1,  where 

(M<J>)(x)  = <P  ' (x)  [Ax  + f Cx,<|) ) ] - B«^(x)  - g (x,4>  (x)  ) . 

Then  as  x - 0,  ||h(x)  - <p  (x)  | | = OC|x|q). 

The  proof  of  Theorem  3 is  the  same  as  that  given  for  the 

finite  dimensional  case  except  that  the  extension  0 : X -*■  Y of  $ 

must  be  defined  so  that  0 (x)  is  in  the  domain  of  B. 

6.3.  Examples 

Example  1.  Consider  the  semilinear  wave  equation 


Vtt  + vt  ' vxx  ' v + £(»  = °>  e CO.11)  x (0 


v = 0 at  x = 0 , Ti 


(6.3.1) 
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where  f is  a C3  function  satisfying  f(v)  ■ v3  ♦ 0(v4)  as 
v -►  0.  We  first  formulate  (6.3.1)  as  an  equation  on  a Hilbert 
space 

Let  Q * (d/dx)'  + 1,  D(Q)  = H2(0,1)  n hJ(0,1).  Then  Q is 
a self-adjoint  operator.  Let  2 = hJ(0,1)  * L2(0,1),  then  (6.3.1) 
can  be  rewritten  as 


w * Cw  ♦ N (w) 


(6.3.2) 


where 


w2 

0 

Cw  * 

Qwrw2 

, N(w)  * 

-f(Wj) 

— — 

_ _ 

Since  C is  the  sum  of  a skew  self-adjoint  operator  and  a bounded 
operator,  C generates  a strongly  continuous  group.  Clearly  N is 
a C^  map  from  2 into  2. 

The  eigenvalues  of  C are  An  * [-1  ± ( 5- 4n* ) 1 ' - ] /2 . 

Aj  = 0 and  all  the  other  eigenvalues  have  real  part  less  than  0. 
The  eigenspace  corresponding  to  the  :ero  eigenvalue  is  spanned  by 
qj  where 


qxU)  * 

To  apply  the  theory  of  Section 
canonical  form.  We  first  note  that 


1 

0 

2 


sin  x. 


, we  must 


put  (6.3.2) 
where 


into 
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q?U)  * 


sin  x 


and  that  all  the  other  eigenspaces  are  spanned  by  elements  of  the 

2 

form  sin  nx,  n > 2 , € IR  . In  particular,  all  other  eigen- 
vectors are  orthogonal  to  and  q2.  Let  X = spanfqj), 

V = span(q1,q2),  Y = span(q2)  @ V1,  then  Z = X ® Y.  The  projection 
P:  Z -*■  X is  given  by 


P = (w1^2)q1 

•2 


(6.3.3) 


where 


w (x)  = | f w. (9 ) sin  e d« . 
J J 0 J 


Let  w = sq^  + y,  s € ]R,  y € Y and  B = (I-P)C,  then  we  can  write 
(6.3.2)  in  the  form 


sqx  = PN(sq1+y) 

y = By  + (I- P)N(sq1+y) 


(6.3.4) 


By  Theorem  1,  (6.3.4)  has  a centre  manifold  y = h(s), 
h(0)  = 0,  h'(0)  = 0,  h:  (-6,6)  Y.  By  Theorem  2,  the  equation 
which  determines  the  asymptotic  behavior  of  solutions  of  (6.3.4) 
is  the  one-dimensional  equation 
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sq j = PN(sq^+h(s) ) . (6.3.5) 

Since  the  nonlinearities  in  (6.3.4)  are  cubic,  h(s)  = 0(s3),  so 
that 

s - '!■  f £(s+0(s3))sin49  d6 
or 

s = j s3  + 0(s4 ) . (6.3.6) 

Hence,  by  Theorem  2,  the  zero  solution  of  (6.3.6)  is  asymptotically 
stable.  Using  the  same  calculations  as  in  Section  1 of  Chapter  3, 
if  s(0)  > 0 then  as  t •+  ® , 

3 ”1/2  -1/2 

s(t)  = (j  t)  + 0(t  i/Z).  (6.3.7) 

Hence,  if  v(x,t)  is  a solution  of  (6.3.1)  with  v(x,0),  vt(x,0) 
small,  then  either  v(x,t)  tends  to  zero  exponentially  fast  or 

v(x,t)  = ±s (t) sin  x + 0(s3)  (6.3.8) 

where  s(t)  is  given  by  (6.3.7). 

Further  terms  in  the  above  asymptotic  expansion  can  be  calculated 

if  we  have  more  information  about  f.  Suppose  that  f(v)  = v3  + 

5 7 

av  + 0(v  ) as  v -*■  0 . In  order  to  calculate  an  approximation  to 
h(s)  set 


■■  mi  n M 


■i 


(M*)U)  - ♦,Cs)PNCsq1H(x))  - BHs)  - ( 1 - P)  N (sq  ♦♦  (s ) ) (6.3.9) 


where  $ : 1R  -*■  Y . To  apply  Theorem  3 we  choose  4>(s)  so  that 
(M*)(s)  - 0(s5).  If  ♦(s)  - 0(s3)  then 


M4>  (s)  - -B$(s)  - (I-P)N(sq1)  ♦ 0(sb) 


-B'f  (s)  - £ s3q2  - 


q ♦ 0(s3) 


(6.3.10) 


where  q(x)  * sin  3x . If 


H«)  - *<\2S'  * 3 9s  3 


then  substituting  (6.3.11)  into  (6.3.10)  we  obtain 


(6.3.11) 


(M$)(s)  ■ aq^s^  ♦ 


2 3 3 3 s3  0 


I,.*,  ' 1 * "*“*■  , 


Hence,  if  a - 3/4,  6 - 1/32,  B,  ■ 0 then  M*(s)  = 0(s5)  so  by 


Theorem  3 


h(s)  - j q2s3  ♦ ^ Q 


qs3  ♦ 0(s5) 


(6.3.12) 


Substituting  (6.3.12)  into  (6.3.5)  we  obtain 


* - ‘t  *3  * <m  * I a)sS  + 0(s?) 


^ ^ • 1 *• 


* , t r t « • 


* i * ? * ^ 

- - *•  v<t  » . . . 


/ / / ✓ 


{ fc  . i u 


xx 


* * J ' ;,J  a,‘‘ 1 k.iv*u  initial  conditions 

° ' * J (*.3.i3j  >s  a model  for  the 

d,(  rod  with  funged  ends,  v being  the 

“ “'»*'*'»*.  The  above  equation  has 

dml  l"  l,H»  titular  he  showed  that  when 
H“  "*  asymptotically  stable. 

' *'  *“  U*'7  "'**  1 "‘^ai  i .'r,J  rquut  ions  have  a zero 

^ i KCUU  4 i tits  -Hot  bo  ( |lu  1 1 o Hi  , . . 

«U  i a\  oi  solutions  depends  on  the 

1111  ">«  '-<•>  "<  .'ecav  of  solutions  of  (6.3.13, 

Hero  so  discuss  the  behavior 

Mt  *‘mU  -olutum-,  ot  s \\)  when  S * i : <Ball. 

v * ‘‘‘  k‘u‘v''‘VW  e sample  we  tovmulate  k{>.  5.151  i5  ir  0rdi-a*-v 

»*  ‘<>1*  ■ for  srsce  derivatives  ari 


^ v'  l It  i d«  V \ V .4 1 v \ <;» 


S V*  t 


‘ ! V V 


■ 1 t< 


: ; 


V . * 


> * 


' ..-e 


s' 
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Then  we  can  write  (6.3.13)  as 

w = Cw  + N(w).  (6.3.14) 

It  is  easy  to  check  that  C generates  a strongly  continuous  group 
on  Z and  that  N is  a C°°  map  from  Z into  Z. 

If  X is  an  eigenvalue  of  C then  we  must  have  a nontrivial 
solution  u(x)  of 


u' ' ' ' - Bu"  + (X+X2)u  = 0 

u(0)  = u" (0)  = u(l)  = u"(l)  = 0. 

An  easy  computation  then  shows  that  X is  an  eigenvalue  of  C if 
and  only  if 


2X  = -1  ± [1  - 4(n4iT4  + Bn2Tt2)]1/2. 

Let  e * tt2B  + it4.  Then  the  eigenvalues  of  C are  X (e)  , where 
2X1(e)  = -1  + [1-4C]1/2,  X2(e)  = -1  - X^e),  and  all  the  rest  of 
the  X (e)  have  real  parts  less  than  zero  for  £ sufficiently 
small.  The  eigenspaces  corresponding  to  Xj(e)  and  X2(e)  are 
spanned  by  q^  and  q2  where 

sin  itx 

while  the  eigenspaces  corresponding  to  Xn(e)  for  n > 2 are 
spanned  by  elements  orthogonal  to  qj  and  q2. 


xx(0 


sin  tix. 


1 

x2(£) 


m 
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Let  X = span Cq j ) , V = spantqj.q,),  Y = span(q2)  © V1,  then 
Z = X © Y and  the  projection  P:  Z X is  given  by 

P 

where 

w . * 2 I w . (6 ) sin  *9  d9 . 

J i o J 

Let  w - sq^  ♦ y,  where  s € 1R  and  y € Y.  Then  we  can  write 
(6.3.14)  in  the  form 

sqi  = Xjfejsqj  + PN(sqj+y) 

? = By  ♦ (I-P)N(sq1+y)  (6.3.15) 

e = 0. 

By  Theorem  1,  (6.3.15)  has  a centre  manifold  y = h(s,e),  |s|  <6, 

< eo ' Using  Ms»e)  = 0(s^+|es|),  if  N2(w)  is  the  second 
component  of  N(w)  then 

N2(sqi*h(s,e))(x)  = — i [f  s ^rc4cos ^rc9d0  Jif^s  sin  fix 

n JO 

♦ 0(s4+|es5|) 

= -s3sin  ttx  * 0(s4+|es3|). 

Hence  PN(sq1+h(s , e) ) = (-s3+0(s4* | eS3 1 ) )qj  , so  that  by  Theorem  2, 


-1,- 


* CAj(c)  - A2(c))-i(w2-X2(e)gi)qi 


n »»fc— 
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the  asymptotic  behavior  of  small  solutions  of  (6.3.15)  is  determined 
by  the  equation 

s = AjCeJs  - s3  + 0(| es3 | + |s4|)  (6.3.16) 

We  can  now  determine  the  asymptotic  behavior  of  small  solutions  of 
(6.3.14).  For  0 < e < <5,  solutions  of  (6.3.14)  are  asymptotically 
stable.  For  -6  < e < 0,  the  unstable  manifold  of  the  origin 
consists  of  two  stable  orbits  connecting  two  fixed  points  to  the 
origin.  (See  Figure  1) 

Problem.  For  e = 0,  show  that  equation  (6.3.16)  can  be  written  as 

3 5 7 

s = -s  - 3s3  + 0 (s  j 


Figure  1 

Phase  Portrait  for  Small  Negative  £ 


■ 


lixample  3.  Consider  the  equation 


utt  + 2u  - u + a2v  + f(u,v)  = 0 
zz  X XX  (6.3.17) 

vtt  + 2vt  - vxx  - u + g(u,v)  = 0 

for  (x,t)  € (0,n)  * (0,°°)  with  u = v = 0 at  x = 0,n,  where 

f(u,v),  g(u,v)  have  a second  order  zero  at  u = v = 0.  For  a = 2, 

we  show  that  the  linearized  problem  has  two  purely  imaginary 

eigenvalues  while  all  the  rest  have  negative  parts.  We  then  use 

centre  manifold  theory  to  reduce  the  prbblem  of  bifurcation  of 

periodic  solutions  to  a two-dimensional  problem. 

. . t 

Let  w = (u , v , u , v)  , then  we  can  write  (6.3  17)  as 


w = Cw  + N(w) 


(6.3.18) 


on  Z = CH2 (0,1)  n Hq (0 , 1 ) ) 2 x (L2 (0 , 1) ) 2 . Let 


t | 

i J 

i (I 


u 

-u  + a v 

XX 

V. 

-v  + u 

* 

XX  J 

1 

H 


If  W is  an  eigenvalue  of  A then  X is  an  eigenvalue  of  C 
2 

where  X + 2X  + p = o and  all  the  eigenvalues  of  C arise  in 
this  way.  An  easy  calculation  shows  that  the  eigenvalues  of  C 
are  given  by 

X * -1  ± (l-n2±ia)1/2,  n = 1,2,.  . . . 
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the  rest  have  negative  real  parts.  Also 


jg  (Re  Xj  (2) ) > 0 

so  that  \^(a)  and  X^(a)  cross  the  imaginary  axis  with  non-zero 
speed.  It  is  now  trivial  to  apply  centre  manifold  theory  to  conclude 
that  for  a - 2 small,  the  behavior  of  small  solutions  is 
determined  by  an  equation  of  the  form 


a -1 
1 a 


s ♦ ,J( s.a) 


(b. 3. 19) 


2 2 
where  s € 1R  , u is  a real  parameter  and  J(s,a)  - 0(s  ) . To 

apply  the  theory  in  Section  2,  Chapter  3,  we  need  to  calculate  the 

quadratic  and  cubic  terms  in  J(s,0).  To  do  this  we  need  to  put 

(6.3.18)  into  canonical  form  and  to  calculate  the  centre  manifold 

when  a = 2.  from  now  on  we  let  a * 2. 

On  the  subspace  (r  sin  nx:  r £ ]R4}  the  operator  C can  be 

represented  by  the  matrix  C where 

n 


0 

0 

n‘ 

1 


0 

0 

-4 

-n‘ 


1 

0 

-2 

0 


0 

1 

0 

„ 2 


Note  that  the  eigenvalues  of  C are  given  by  the  eigenvalues  of 


for  n *=  1,2,.. 


To  put  lb. 3. 18)  into  canonical  form  we 


first  find  a basis  which  puts  Cj  into  canonical  form.  Calculations 


i 


F 

I 

I 


\ 


i» 

i 


l> 


show  that  if 


-2 

0 

0 

-2 

0 

1 

-1 

0 

0 

<2  = 

-2 

«3  = 

-2 

<*4  = 

4 

-1 

0 

2 

1 

then  C2^1  = ’q2’  Ciq”>  = qi»  Ciq3  = '2^3  + q4  » 
Let  w = Qz  where  Q = [q x ,q2  ,q 3 ,q4  ] , then  we 
as 

2 = Q 1CQz  + Q"1N(Qw). 

Let  X = { [s^s^O,  0 I Tv  : Sl,s2  € IR } , V = {(0,0 
Y = V @ [X  0 V]"*"  where  V(x)  = sin  x.  Then  Z 
projection  P on  X along  Y is  given  by 


Pw 


Wj  (x) 


[ w. (0)sin  OdO . 
JO  J 


Let 

(6.3.20) 


z = (s1,s2,0,0]Tv 
in  the  form 


♦ y, 


st  6 JR,  y € Y, 


Clq4  = ’q3  ‘ 2q4' 
can  rewrite  (6.3.18) 


(6.3.20) 

’ r 1 » r 2 ] T|f  : rltr2  € 3R } , 
= X © Y and  the 


then  we  can  write 
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0 -1 

1 0 


s ♦ <PQ'1N(Qw)  ,♦> 


y = By  + (I *P)Q*  xN(Qw) 


(6.3.21) 


where  B * (I-PJQ^CQ,  * * [1,1,0,0]T¥  and  s * [s^s^1.  By 
Theorem  1 (6.3.21)  has  a centre  manifold  y = h(s). 

From  w = Qz  we  have  that  Wj  = -2Zj  - 2z4  and  w2  =■  z,  - z 
Let  Q*1  = (tij), 


F(z) 

K 

tl 3 *14 

f(*2Zj  - 2z4  , z2  * z3) 

G(z) 

i23  l2£ 

8(*2Zj  - 2z4 , z2  - Zj) 

Suppose  that  F(z)  = F3(z)  ♦ 0(|z|4)  and  G(z)  « G3(z)  ♦ 0(|z|4) 
where  F3  and  F4  are  homogeneous  cubics.  Then  if  J^s),  J,(s) 
denote  the  first  two  components  of  <PQ~ lN(Qw) ,<P>  on  the  centre 
manifold , 


2 f71 

Jl(s^  = f J FjUjSin  6,  s2sin  0,O,O)d0  ♦ 0(|s|4) 


with  a similar  expression  for  J2(s).  Hence,  on  the  centre  manifold 


0 -1 

S1 

♦ 

Jj(s) 

l 0 

s: 

J2(s) 

_ _ 

_ — 

L J 

(6.3.22) 


and  we  can  apply  the  theory  given  in  Section  2 of  Chapter  3.  If  the 
constant  K associated  with  (6.3.22)  is  zero  (see  Section  2 of 


"" 
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Chapter  3 tor  the  definition  of  K)  then  the  above  procedure  gives 
no  information  and  we  have  to  calculate  higher  order  terms. 

If  F(z)  = F2U)  + Fj(z)  ♦ OC  | z | 4 ) and  G(z)  = G2(z)  ♦ 

G3(z)  ♦ 0 C I z I 4 ) where  and  G^  are  homogeneous  quadratics 

then  the  calculation  of  the  nonlinear  terms  is  much  more  complicated. 
On  the  centre  manifold  Zj  = SjY  ♦ 0(s2),  z2  = + 0(s2), 

2 2 7 

2 3 = 0(s"),  and  z4  * Ots*").  The  terms  of  order  s^  make  a contri- 
bution to  the  cubic  terms  in  PQ_1NCQw).  Hence,  we  need  to  find  a 
quadratic  approximation  to  h(s).  This  is  straightforward  but  rather 
complicated  so  we  omit  the  details. 

F:xample  4.  Consider  the  equations 

ut  = Duxx  + * a2v  + 2Auv  + u 2 v ♦ BA'*u2 

Vt  = 6Dvxx  ' Bu  ' a2v  * 2Auv  - u2v  - BA* Ju2  (6.3.23) 

u=v=0  at  x * 0,1, 

where  A,B,0,D  are  positive.  The  above  equations  come  from  a 
simplified  model  of  a chemical  reaction  with  u + A and  v + BA’^ 
as  the  chemical  concentrations  [3,10]. 

We  study  (6.3.23)  on  Z - (hJ(0,tt))2. 


Set 


1 
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w 


D — *-  ♦ ( B • 1 ) 
dxz 


-B 


0D  -i-T  - A2 


dx 


N (w) 


(2Auv  u2v  ♦ BA'1!!2) 


then  we  can  write  (6.3.23)  as 


w * Cw  ♦ N(w) 


(6.3.24) 


We  analyze  the  situation  in  which  for  some  value  of  the  para- 
meters A,B,9,D,  C has  two  zero  eigenvalues  such  that  the  restriction 
of  C to  the  zero  eigenspace  has  a Jordan  block.  The  bifurcation 
of  static  solutions  when  C has  two  zero  eigenvalues  and  the 
restriction  of  C to  the  zero  eigenspace  is  zero,  has  been  studied 
in  [21,22]. 


On  the  subspace  (r  sin  n^x,  r € ]RZ)  the  operator  C can  be 


. I 


<4 


I 
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represented  by  the  matrix 


n 


-n2Ti2D  ♦ B - 1 


-B 


-t)Dn2"2  - A2 


The  eigenvalues  of  C are  given  by  the  eigenvalues  of  Cn  for 


n = 1,2 We  suppose  that  two  of  the  eigenvalues  of  C are 

negative  while  all  the  rest  have  negative  real  parts.  For 


simplicity  we  assume  that  the  eigenvalues  of  are  zero. 


If  Cj  is  to  have  two  zero  eigenvalues  then 


trace  (Cj ) = B - 1 - u2n  - A2  - Qtt2d  = 0 


det  (Cj ) = A2B  - (B-l-n2D)  (A2+0"2D)  ■=  0, 


(6.3.25) 


1 

- 


1 


1 


r 

' 


We  make  the  following  hypotheses: 

(HI)  There  exists  A0*Bo’eo,D0  such  that  Cb.3.25)  is  satisfied 
and  the  real  parts  of  the  rest  of  the  eigenvalues  of  C 
are  negative. 

(H2)  For  A,B,0,D  in  a neighborhood  of  Ag  >Bg  »%  »Dp  * we  can 
parametrize  trace(C1)  and  det(Cj)  by 


trace (C x ) = e2>  det(Cx)  * -Cj.  (6.3.26) 
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The  first  hypothesis  is  satisfied,  for  example,  if  DQ  > 0 


and 


% ■ »•  Ao  • "V  Bo  ■ »*”2v2 


(6.3.27) 


I f we  vary  B and  8 and  keep  A = Ag , D = Dg , then  the  mapping 
(B,6)  -*•  (det(C^),  trace(C^))  has  a non-zero  Jacobian  at  B = BQ , 
0 = if  Ao»Bo,0O  are  8*ven  by  (6.3.27),  so  by  the  implicit 
function  theorem,  (H2)  is  satisfied.  In  order  to  simplify 
calculations  we  assume  (6.3.27)  from  now  on. 

Let  X = (sV:  s £IR^},  Y = where  'j'(x)  = sin  *x,  then 
Z = X 0 Y.  By  Theorem  1,  the  system 


i 


w = Cw  ♦ N(w) 

e = o 


has  a centre  manifold  h:  (neighborhood  of  X * IR^)  -►  Y,  where  we 
have  written  e * (e1.e2^*  0n  t*ie  centre  manifold,  the  equation 
reduces  to 


s * (^(ejs  ♦ N(s<l'+h(s , e) ) 


(6.3.28) 


where  s = [Sj,s2]  and 


N.(z)  * 2 f N.(z(0))sin  ir6d0 
1 Jn1 


We  treat  the  linear  and  nonlinear  parts  of  (6.3.28)  separately, 


n>  tm 1 — i 


~“1 
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If  q!  = [lfPlT,  q2  = [ - P » 1 J 1 Where  p 
Cj(0)qj  = 0 and  C^COJq^  = (A"+B)qj.  Let  Q 

0 A2  + B 

0 0 

L- 

Let  Q^CjCOQ  - T - and  let 


■ A’ 1 B1^2 
A0  B0  * 

[q1»q2l » 


Q Cj(0)Q  = 


M(e)  = 


1 

r 

o 

f* 

1 — 
o 

hi 

t12 

1 o 

a2+b 

0(e) 


Then  for  e sufficiently  small,  M(e)  is  nonsingular  and 


M(e)Q'1C1(e)QM'1(e)  . 


0 1 
■det(T)  trace(T) 


i 

0 

1 

0 

1 

-detCCj) 

trace (C^) 

£1 

£2 

by  (b.3.26).  Let 

s » QM'1 

(e)r,  then  (6.3. 

,28) 

becom 

0 1 

ei  £; 


r ♦ M(e)Q'1N(QM‘1(e)r<|»  ♦ h(s,e)) 


then 

then 


r « 


(t).  3. 29) 
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To  check  the  hypotheses  of  Section  9 of  Chapter  4 we  only 

need  calculate  the  nonlinear  terms  when  e = 0.  Using  the  fact 

2 

that  h(s,0)  * 0(s  ),  routine  calculations  show  that 


M(0)Q*  1R(QM" 1 (O)r'i'  ♦ h(s  ,0) ) = [px ,P2J lR(rx ,r2)  * 0(|rj|  ♦ |r2|J) 


where 


pi  = (^+p2 *) " 1 (i *p)  » p2  » 3I  U+pW+PMA^B) 


R(r1,r2)  = + 01 2 r 1 r 2 + a3r2 


ri  . 91,1/2 


a = BA  A - 2Bi/<',  0 = 2 (A2+B)  (B3//2A'2+A-BA'1) 


Using  (6.3.27),  = (tt2Dq)  *(l-n4 *Dp)  so  that  is  non-zero 

2 2 

1 f 11  Dp  ♦ 1 . We  assume  that  * 1 from  now  on.  Note,  a lso 

2 -I 

that  since  1 ♦ p = - (n  Dp)  , we  have  that  p2  is  non-zero. 

To  reduce  (6.3.29)  to  the  form  given  in  Section  9,  Chapter  4 

we  make  the  substitution 


-1, 


P * (I-p^  A(£))r 


(6.3.30) 


Substituting  (6.3.30)  into  (6.3.29)  and  using  the  above  calculations 
we  obtain 


0 1 


£1  £2 


P + F(P,e) 


(6.3.31) 
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32F,(0,0) 

V2>  6 ’ -3ppts 2Vl  * *2*2 

0. 

ecked  that  a.  * 0 . For  most  values  of 
only  zero  when  Dq  is  a solution  of  a 
certain  algebraic  equation].  If  B * 0 then  we  can  apply  the 
theory  given  in  Section  9,  Chapter  4 to  obtain  the  bifurcation  set 
for  (6.3.23).  If  B = 0 then  the  theory  given  in  Section  9, 
Chapter  4 still  gives  us  part  of  the  bifurcation  set;  the  full 
bifurcation  sec  would  depend  on  higher  order  terms. 

Remark . If  we  vary  9 in  (6.3.23)  the  theory  given  in  Section  l 
does  not  apply  since  the  map  (9  ,v)  0vxx  is  not  even  defined  on 
the  whole  space.  However,  it  is  easy  to  modify  the  results  of 
Section  1 to  accommodate  the  above  situation.  (See,  for  example, 
Exercises  1-2  in  Section  3.4  of  [30].) 


9“F7(0,0) 


9p: 


3Fj(0,0) 

2 

9p: 


We  have  already  ch 
D0,B  is  nonzero  IB  is 
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